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BIOINFORMATICALLY DETECTABLE GROUP OF NOVEL 
REGULATORY GENES AND USES THEREOF 

BACKGROUND OF THE INVENTION 

FIELD OF THE INVENTION 

The present invention relates to a group of bioinfonnatically detectable novel 
genes, here identified as Genomic Address Messenger or GAM genes, which are believed to be 
related to the micro RNA (miRNA) group of genes. 



DESCRIPTION OF PRIOR ART 

MIR genes are regulatory genes encoding MicroRNA's (miRNA), short ~22nt 
non-coding RNA's, found in a wide range of species, believed to function as specific gene 
translation repressors, sometimes involved in cell-differentiation. Some 110 human MIR genes 
have been detected by laboratory means. Over the past 6 months, the need for computerized 
detection of MIR genes has been recognized, and several informatic detection engines have been 
reported (Lim, 2 003; Grad, 2 003; Lai, 2 003). Collectively these informatic d etection engines 
found 38 more human MIR genes which were later confirmed in zebrafish, 14 human MIRs 
which were confirmed in human, and 55 postulated human MIRs which could not be confirmed 
by laboratory (Lim, 2003). Extensive efforts to identify novel MIR genes using c onventional 
biological detection techniques such as massive cloning and sequencing efforts, and several 
bioinformatic detection attempts, have led leading researchers in the field to the conclusion that 
the total number of human MIR genes is between 200 to 255 (Lau, 2003; Lim 2003 Science; 
Lim, 2003 Genes Dev). Recent studies postulate that the number of MIR genes may be higher 
(Grad, 2003; Krichevsky, 2003). 

The ability to detect novel MIR genes is limited by the methodologies used to 
detect such genes. All MIR genes identified so far either present a visibly discernable whole 
body phenotype, as do Lin-4 and Let-7 (Wightman ,1993; Reinhart , 2000), or produce sufficient 
quantities of RNA so as to be detected by the standard molecular biological techniques. 

Initial studies reporting MIR genes (Bartel, 2001; Tuschl, 2001) discovered 93 
MIR genes in several species, by sequencing a limited number of clones (300 by Bartel and 100 
by Tuschl) of small segments (i.e. size fractionated) RNA. MiRNA encoded by MIR genes 
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detected in these studies therefore, represent the more prevalent among the miRNA gene family, 
and can not be much rarer than 1% of all small ~20nt-long RNA segments. 

Current methodology has therefore been unable to detect micro RNA genes (MIR 
genes) which either do not present a visually discernable whole body phenotype, or are rare (e.g. 
rarer than 0.1% of all size fractionated ~20nt-long RNA segments expressed in the tissues 
examined), and therefore do not produce significant enough quantities of RNA so as to be 
detected by standard biological techniques. 



BRIEF DESCRIPTION OF SEQUENCE LISTING, LARGE TABLES AND COMPUTER 

PROGRAM LISTING 

Sequence listing, large tables related to sequence listing, and computer program 
listing are filed under section 8 01 (a)(i)onan electronic medium in computer readable form, 
attached to the present invention, and are hereby incorporated by reference. Said sequence 
listing, large tables related to sequence listing, and computer program listing are submitted on a 
CD-ROM (Operating system: MS-Windows), entitled SEQUENCE LISTING AND LARGE 
TABLES, containing files the names and sizes of which are as follows: 

Sequence listing comprising 548,156 genomic sequences, is filed under section 
801 (a)(i) on an electronic medium in computer readable form, attached to the present invention, 
and is hereby incorporated by reference. Said sequence listing is contained in a self extracting 
compressed file named SEQ_LIST.EXE (9,813KB). Compressed file contains 1 file named 
SEQ_LIST.TXT (82,859KB). 

Large t ables r elating t o genomic s equences are stored i n 7 s elf e xtracting files, 
each comprising a respective one of the following table files: TABLE1.TXT (839KB); 
TABLE2.TXT (4348KB); TABLE3.TXT (715KB); TABLE4.TXT (545728KB); TABLE5.TXT 
(1523981KB); TABLE6.TXT (373091KB); and TABLE7.TXT (125KB). 

Computer program listing of a computer program constructed and operative in 
accordance with a preferred embodiment of the present invention is enclosed on an electronic 
medium in computer readable form, and is hereby incorporated by reference. The computer 
program listing is contained in a self extracting compressed file named COMPUTER 
PROGRAM LISTDSTG.EXE (168KB). Compressed file contains 5 files, the name and sizes of 
which are as follows: AUXILARYJFILES.TXT (117K); BINDING_SITE_SCORING.TXT 
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(17K); EDITJDISTANCE.TXT (144K); FIRST-K.TXT (96K); HAIRPIN_PREDICTIO.TXT 
(47K); and TWO_PHASED_PREDICTOR.TXT (74K). 
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SUMMARY OF THE INVENTION 

The present invention relates to a novel group of regulatory, non-protein coding 
genes, which are functional in specifically inhibiting translation of target proteins. Each gene in 
this novel group of genes, here identified as GAM or Genomic Address Messengers, specifically 
inhibits translation of one of more other 'target' genes by means of complementary hybridization 
of a segment of the RNA transcript encoded by GAM, to an inhibitor site located in an 
untranslated region (UTR) of the mRNA of the one or more target' genes. 

In various preferred embodiments, the present invention seeks to provide 
improved method and system for specific modulation of expression of specific known 'target' 
genes involved in significant human diseases, and improved method and system for detection of 
expression of novel genes of the present invention, which modulate these target genes. 

Accordingly, the invention provides several substantially pure DNAs (e.g., 
genomic DNA, cDNA or synthetic DNA) each encoding a novel gene of the GAM group of gene, 
vectors comprising the DNAs, probes comprising the DNAs, a method and system for selectively 
modulating translation of known 'target' genes utilizing the vectors, and a method and system for 
detecting expression of known 'target' genes utilizing the probe. 

By "substantially pure DNA" is meant DNA that is free of the genes which, in the 
naturally-occurring genome of the organism from which the DNA of the invention is derived, 
flank the genes discovered and isolated by the present invention. The term therefore includes, for 
example, a recombinant DNA which is incorporated into a vector, into an autonomously 
replicating plasmid or virus, or into the genomic DNA of a prokaryote or eukaryote at a site other 
than its natural site; or which exists as a separate molecule (e.g., a cDNA or a genomic or cDNA 
fragment produced by PCR or restriction endonuclease digestion) independent of other 
sequences. It also includes a recombinant DNA which is part of a hybrid gene encoding 
additional polypeptide sequence. 

"Inhibiting translation" is defined as the ability to prevent synthesis of a specific 
protein encoded by a respective gene, by means of inhibiting the translation of the mRNA of this 
gene. "Translation inhibitor site" is defined as the minimal DNA sequence sufficient to inhibit 
translation. 
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There is thus provided in accordance with a preferred embodiment of the present 
invention a b ioinformatically d etectable n ovel g ene e ncoding s ubstantially pureD NA w herein 
RNA encoded by the bioinformatically detectable novel gene is about 18 to about 24 nucleotides 
in length, and originates from an RNA precursor, which RNA precursor is about 50 to about 120 
nucleotides in length, a nucleotide sequence of a first half of the RNA precursor is a partial or 
accurate inversed reversed sequence of a nucleotide sequence of a second half thereof, a 
nucleotide sequence of the RNA encoded by the novel gene is a partial or accurate inversed 
reversed sequence of a nucleotide sequence of a binding site associated with at least one target 
gene, the novel gene cannot be detected by either of the following: a visually discernable whole 
body phenotype, and detection of 99.9% of RNA species shorter than 25 nucleotides expressed in 
a tissue sample, and a function of the novel gene is bioinformatically deducible. 

There is still further provided in accordance with another preferred embodiment of 
the present invention a bioinformatically detectable novel gene encoding substantially pure DNA 
wherein RNA encoded by the bioninformatically d etectable novel gene includes a plurality of 
RNA sections, each of the RNA sections being about 50 to about 120 nucleotides in length, and 
including an RNA segment, which RNA segment is about 18 to about 24 nucleotides in length, a 
nucleotide sequence of a first half of each of the RNA sections encoded by the novel gene is a 
partial or accurate inversed reversed sequence of nucleotide sequence of a second half thereof, a 
nucleotide sequence of each of the RNA segments encoded by the novel gene is a partial or 
accurate inversed reversed sequence of the nucleotide sequence of a binding site associated with 
at least one target gene, and a function of the novel gene is bioinformatically deducible from the 
following data elements: the nucleotide sequence of the RNA encoded by the novel gene, a 
nucleotide sequence of the at least one target gene, and function of the at least one target gene. 

There is additionally provided in accordance with another preferred embodiment 
of the p resent i nvention a b ioinformatically d etectable n ovel gene e ncoding s ubstantially p ure 
DNA wherein RNA encoded by the bioinformatically detectable novel gene is about 18 to about 
24 nucleotides in length, and originates from an RNA precursor, which RNA precursor is about 
50 to about 120 nucleotides in length, a nucleotide sequence of a first half of the RNA precursor 
is a partial or accurate inversed reversed sequence of a nucleotide sequence of a second half 
thereof, a nucleotide sequence of the RNA encoded by the novel gene is a partial or accurate 
inversed reversed sequence of a nucleotide sequence of a binding site associated with at least one 
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target gene, a function of the novel gene is modulation of expression of the at least one target 
gene, and the at least one target gene does not encode a protein. 

There is moreover provided in accordance with another preferred embodiment of 
the present invention a bioinformatically detectable novel gene encoding substantially pure DNA 
wherein the bioinformatically detectable novel gene does not encode a protein, RNA encoded by 
the bioinformatically detectable novel gene is maternally transferred by a cell to at least one 
daughter cell of the cell, a function of the novel gene includes modulation of a cell type of the 
daughter cell, and the modulation is bioinformatically deducible. 

There is further provided in accordance with another preferred embodiment of the 
present invention a bioinformatically detectable novel gene encoding substantially pure DNA 
wherein the bioinformatically detectable novel gene does not encode a protein, a function of the 
novel gene is promotion of expression of the at lease one target gene, and the at least one target 
gene is bioinformatically deducible. 

Still further in accordance with a preferred embodiment of the present invention 
the function of the novel gene is bioinformatically deducible from the following data elements: 
the nucleotide sequence of the RNA encoded by the bioinformatically detectable novel gene, a 
nucleotide sequence of the at least one target gene, and a function of the at least one target gene. 

Additionally in accordance with a preferred embodiment of the present invention 
the RNA encoded by the novel gene complementarity binds the binding site associated with the 
at least one target gene, thereby modulating expression of the at least one target gene. 

Moreover in accordance with a preferred embodiment of the present invention the 
binding site associated with at least one target gene is located in an untranslated region of RNA 
encoded by the at least one target gene. 

Further in accordance with a preferred embodiment of the present invention the 
function of the novel gene is selective inhibition of translation of the at least one target gene, 
which selective inhibition includes complementary hybridization of the RNA encoded by the 
novel gene to the binding site. 

Still further in accordance with a preferred embodiment of the present invention, 
the present invention includes a vector including the DNA 
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Additionally in accordance with a preferred embodiment of the present invention, 
the present invention includes a method of selectively inhibiting translation of at least one gene, 
including introducing the vector into a cell. 

Moreover in accordance with a preferred embodiment of the present invention the 
introducing includes utilizing RNAi pathway. 

Further in accordance with a preferred embodiment of the present invention, the 
present invention includes a gene expression inhibition system including: the vector and a vector 
inserter, functional to insert the vector of claim 10 into a cell, thereby selectively inhibiting 
translation of at least one gene. 

Still further in accordance with a preferred embodiment of the present invention, 
the present invention includes a probe including the DNA 

Additionally in accordance with a preferred embodiment of the present invention, 
the present invention includes a method of selectively detecting expression of at least one gene, 
including using the probe 

Moreover in accordance with a preferred embodiment of the present invention, the 
present invention includes a gene expression detection system including: the probe, and a gene 
expression detector functional to selectively detect expression of at least one gene. 

BRIEF DESCRIPTION OF DRAWINGS 

Fig. 1 is a simplified diagram illustrating the genomic differentiation enigma that 
the present invention addresses; 

Figs. 2,3 and 4 are schematic diagrams which when taken together provide an 
analogy that illustrates a conceptual model of the present invention, addressing the genomic 
differentiation enigma; 

Figs. 5A and 5B are schematic diagrams, which when taken together illustrate a 
'genomic records" concept of the conceptual model of the present invention, addressing the 
genomic differentiation enigma; 
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Fig. 6 is a schematic diagram illustrating a "genomically programmed cell 
differentiation" concept of the conceptual model of the present invention, addressing the genomic 
differentiation enigma; 

Fig. 7 is a schematic diagram illustrating a "genomically programmed cell-specific 
protein expression modulation concept of the conceptual model of the present invention, 
addressing the genomic differentiation enigma; 

Fig. 8 is a simplified diagram illustrating a mode by which genes of a novel group 
of genes of the present invention, modulate expression of known target genes; 

Fig. 9 is a simplified block diagram illustrating a bioinformatic gene detection 
system capable- of detecting genes of the novel group of genes of the present invention, which 
system is constructed and operative in accordance with a preferred embodiment of the present 
invention; 

Fig. 10 is a simplified flowchart illustrating operation of a mechanism for training 
of a computer system to recognize the novel genes of the present invention, which mechanism is 
constructed and operative in accordance with a preferred embodiment of the present invention; 

Fig. 1 1 A is a simplified block diagram of a non-coding genomic sequence detector 
constructed and operative in accordance with a preferred embodiment of the present invention; 

Fig. 1 IB is a simplified flowchart illustrating operation of a non-coding genomic 
sequence detector constructed and operative in accordance with a preferred embodiment of the 
present invention; 

Fig. 12 A is a simplified block diagram of a hairpin detector constructed and 
operative in accordance with a preferred embodiment of the present invention; 

Fig. 12B is a simplified flowchart illustrating operation of a hairpin detector 
constructed and operative in accordance with a preferred embodiment of the present invention; 

Fig. 13A is a simplified block diagram of a dicer-cut location detector constructed 
and operative in accordance with a preferred embodiment of the present invention; 

Fig. 13B is a simplified flowchart illustrating training of a dicer-cut location 
detector constructed and operative in accordance with a preferred embodiment of the present 
invention; 
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Fig. 13C is a simplified flowchart illustrating operation of a dicer-cut location 
detector constructed and operative in accordance with a preferred embodiment of the present 
invention; 

Fig. 14A is a simplified block diagram of a target-gene binding-site detector 
constructed and operative in accordance with a preferred embodiment of the present invention; 

Fig. 14B is a simplified flowchart illustrating operation of a target-gene binding- 
site detector constructed and operative in accordance with a preferred embodiment of the present 
invention; 

Fig. 15 is a simplified flowchart illustrating operation of a function & utility 
analyzer constructed and operative in accordance with a preferred embodiment of the present 
invention; 

Fig. 16 is a simplified diagram describing a novel bioinformatically detected 
group of regulatory genes, referred to here as Genomic Record (GR) genes, each of which 
encodes an ^operon-like* cluster of novel miRNA-like genes, which in turn modulate expression 
of one or more target genes; 

Fig. 17 is a simplified diagram illustrating a mode by which genes of a novel 
group of operon-like genes of the present invention, modulate expression of other such genes, in 
a cascading manner; 

Fig. 18 is a block diagram illustrating an overview of a methodology for finding 
novel genes and novel operon-like genes of the present invention, and their respective functions; 

Fig. 19 is a block diagram illustrating different utilities of novel genes and novel 
operon-like genes, both of the present invention; 

Figs. 20A and 20B are simplified diagrams, which when taken together illustrate a 
mode of gene therapy applicable to novel genes of the present invention; 

Fig. 21 is a table summarizing laboratory validation results which validate efficacy 
of a bioinformatic gene detection system constructed and operative in accordance with a 
preferred embodiment of the present invention; 

Fig. 22 is a picture of laboratory results validating the expression of 25 novel 
genes detected by a bioinformatic gene detection engine constructed and operative in accordance 
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with a preferred embodiment of the present invention, thereby validating the efficacy of the gene 
detection engine of the present invention; 

Fig. 23A is a schematic representation of an c operon like' cluster of novel gene 
hairpin sequences detected bioinformatically by a bioinformatic gene detection engine 
constructed and operative in accordance with a preferred embodiment of the present invention, 
and non-GAM hairpin useful as negative controls thereto; 

Fig. 23B is a schematic representation of secondary folding of hairpins of the 
operon-like cluster of Fig. 23 A; 

Fig. 23 C is a picture of laboratory results demonstrating expression of novel genes 
of Figs. 23 A and 23B, and lack of expression of the negative controls, thereby validating efficacy 
of bioinformatic detection of GAM genes and GR genes of the present invention, by a 
bioinformatic gene detection engine constructed and operative in accordance with a preferred 
embodiment of the present invention; 

Fig. 24A is an annotated sequence of EST72223 comprising known miRNA gene 
MIR98 and novel gene GAM25, both detected by the gene detection system of the present 
invention; and 

Figs. 24B and 24C are pictures of laboratory results demonstrating laboratory 
confirmation of expression of known gene MIR98 and of novel bioinformatically detected gene 
GAM25 respectively, both of Fig. 24A, thus validating the bioinformatic gene detection system 
of the present invention. 

DETAILED DESCRIPTION OF DRAWINGS 

Reference is now made to Fig. 1 which is a simplified diagram providing a 
conceptual explanation of a genomic differentiation enigma, which the present invention 
addresses. 

Fig. 1 depicts different cell types in an organism, such as CARTILAGE CELL, 
LIVER CELL, FIBROBLAST CELL and BONE CELL all containing identical DNA, and 
deriving from the initial FERTILIZED EGG CELL, and yet each of these cells expressing 
different proteins, and hence acquiring different shape and function. 
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The present invention proposes that the inevitable conclusion from this constraint 
is, however, strikingly simple: the coding system used must be modular . It must comprise 
multiple modules, or records, one for each cell-type, and a mechanism whereby each cell at its 
inception is instructed which record to open, and behaves according to instructions in that record. 

This modular code concept is somewhat difficult to grasp, since we are strongly 
habituated to viewing things from an external viewpoint. An architect, for example, looks at a 
blueprint of a building, which details exactly where each element (block, window, door, electrical 
switch, etc.) is to be placed relative to all other elements, and then instructs builders to place 
these elements in their designated places. This is an external viewpoint: the architect is external 
to the blueprint, which itself is external to the physical building, and its different elements. The 
architect may therefore act as an 'external organizing agent': seeing the full picture and the 
relationships between all elements, and being able to instruct from the outside where to place 
each of them. 

Genomics differentiation coding evidently works differently, without any such 
external organizing agent: It comprises only one smart block (the first cell), which is the 
architect and the blueprint, and which continuously duplicates itself, somehow knowing when to 
manifest itself as a block and when as a window, door, or electrical switch. 

Reference is now made to Figs. 2 through 4 which are schematic diagrams which 
when taken together provide an analogy that illustrates a conceptual model of the present 
invention, addressing the genomic differentiation enigma. 

Reference is now made to Fig. 2A. Imagine a very talented chef, capable of 
preparing any meal provided he is given specific written cooking instructions. This chef is 
equipped with two items: (a) a thick recipe book, and (b) a small note with a number scribbled on 
it. The book comprises multiple pages, each page detailing how to prepare a specific meal. The 
small note indicates the page to be opened, and therefore the meal to be prepared. The chef looks 
at the page-number written on the note, opens the recipe book at the appropriate page, and 
prepares the meal according to the written instructions on this page. As an example, Fig. 2A 
depicts a CHEF holding a note with the number 12 written on it, he opens the book on page 12, 
and since that page contains the recipe for preparing BREAD, the CHEF prepares a loaf of 
BREAD. 



11 



WO 2004/044123 PCT/IL2003/000970 

Reference is now made to Fig. 2B, which depicts two identical chefs, CHEF A 
and CHEF B, holding an identical recipe book. Despite their identity, and the identity of their 
recipe book, since CHEF A holds a note numbered 12, and therefore opens the book on page 12 
and prepares BREAD, whereas CHEF B holds a note numbered 34 and therefore opens the book 
on page 34 and prepares a PIE. 

Reference is now made to Fig. 3. Imagine the chef of the analogy is also capable 
of duplicating himself once he has finished preparing the specified meal. The format of the book 
is such that at the bottom of each page, two numbers are written. When he has finished preparing 
the meal specified on that page, the chef is trained to do the following: (i) divide himself into two 
identical duplicate chefs, (ii) duplicate the recipe book and hand a copy to each of his duplicate 
chefs, and (iii) write down the two numbers found at the bottom of the page of the meal he 
prepared, on two small notes, handing one note to each of his two duplicate chefs. 

Each of the two resulting duplicate chefs are now equipped with the same book, 
and have the same talent to prepare any meal, but since each of them received a different note, 
they will now prepare different meals. 

Fig. 3 depicts CHEF A holding a recipe book and receiving a note numbered 12. 
CHEF A therefore opens the book on page 12 and prepares BREAD. When he is finished 
making bread, CHEF A performs the following actions: (i) divides himself into two duplicate 
chefs, designated CHEF B and CHEF C, (ii) duplicates his recipe book handing a copy to each of 
CHEF B and CHEF C, (iii) writes down the numbers found at the bottom of page 12, numbers 34 
and 57, on two notes, handing note numbered 34 to CHEF B and note numbered 57 to CHEF C. 

Accordingly, CHEF B receives a note numbered 34 and therefore opens the recipe 
book on page 34 and prepares PIE, whereas CHEF C receives a note numbered 57 and therefore 
opens the book on page 57 and therefore prepares RICE. 

It is appreciated that while CHEF A, CHEF B & CHEF C are identical and hold 
identical recipe books, they each prepare a different meal. It is also appreciated that the meals 
prepared by CHEF B and CHEF C are determined CHEF A, and are mediated by the differently 
numbered notes passed on from CHEF A to CHEF B and CHEF C. 

It is further appreciated that the mechanism illustrated by Fig. 3 enables an 
unlimited lineage of chefs to divide into duplicate, identical chefs and to determine the meals 
those duplicate chefs would prepare. For example, having been directed to page 34, when CHEF 
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B divides into duplicate chefs (not shown), he will instruct its two duplicate chefs to prepare 
meals specified on pages 14 and 93 respectively, according to the numbers at the bottom of page 
34 to which he was directed. Similarly, CHEF C will instruct its duplicate chefs to prepare meals 
specified on pages 21 and 46 respectively, etc. 

Reference is now made to Fig. 4. Imagine that the cooking instructions on each 
page of the recipe book are written in shorthand format: The main meal-page to which the chef 
was directed by the scribbled note, merely contains a list of numbers which direct him to multiple 
successive pages, each specifying how to prepare an ingredient of that meal. 

As an example, Fig. 4 depicts CHEF A of FIGS 2 and 3, holding a recipe book 
and a note numbered 12. Accordingly, CHEF A opens the recipe book on page 12, which details 
the instructions for preparing BREAD. However, the 'instructions' on making BREAD found on 
page 12 comprise only of 3 numbers, 18, 7 and 83, which 'refer' CHEF A to pages detailing 
preparation of the ingredients of BREAD - FLOUR, MILK and SALT, respectively. 

As illustrated in Fig. 4, turning from the main 'meal page' ( e.g. 12) to respective 
'ingredients pages' (e.g. pages 18, 7 & 83) is mediated by scribbled notes with the page-numbers 
written on them. In this analogy, the scribbled notes are required for seeking the target pages to 
be turned to - both when turning to main 'meal pages' (e.g. page 12), as well as when turning to 
'ingredient pages' (e.g. pages 18, 7 & 83). 

The chef in the given analogy, schematically depicted in FIGS 2 through 4, 
represents a cell; the thick recipe book represents the DNA; preparing a meal in the given 
analogy represents the cell manifesting itself as a specific cell-type; and ingredients of a meal 
represent proteins expressed by that cell-type. Like the chef equipped with the thick recipe book 
in the given analogy, all cells in an organism contain the same DNA and are therefore each 
potentially capable of manifesting itself as any cell-type, expressing proteins typical of that cell 
type. 

Reference is now made to Figs. 5A and 5B which are schematic diagrams, which 
when taken together illustrate a 'genomic records' concept of the conceptual model of the 
present invention, addressing the genomic differentiation enigma. 

The Genomic Records concept asserts that the DNA (the thick recipe book in the 
illustration) comprises a very large number of Genomic Records (analogous to pages in the recipe 
book), each containing the instructions for differentiation of a different cell-type, or 
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developmental process. Each Genomic Record is headed by a very short genomic sequence 
which functions as a 'Genomic Address' of that Genomic Record (analogous to the page number 
in the recipe book). At its inception, in addition to the DNA, each cell also receives a short RNA 
segment (the scribbled note in the illustration). This short RNA segment binds complementarity 
to a 'Genomic Address' sequence of one of the Genomic Records, thereby activating that 
Genomic Record, and accordingly determining the cell's-fate (analogous to opening the book on 
the page corresponding to the number on the scribbled note, thereby determining the meal to be 
prepared). 

Reference is now made to Fig. 5A. Fig 5A illustrates a CELL which comprises a 
GENOME. The GENOME comprises a plurality of GENOMIC RECORDS, each of which 
correlates to a specific cell type (for clarity only 6 sample genomic records are shown). Each 
genomic record comprises genomic instructions on differentiation into a specific cell-type, as 
further elaborated below with reference to Fig. 7. At cell inception, the CELL receives a 
maternal short RNA segment, which activates one of the GENOMIC RECORDS, causing the cell 
to differentiate according to the instructions comprised in that genomic record. As an example, 
Fig. 5A illustrates reception of a maternal short RNA segment designated A 5 and outlined by a 
broken line, which activates the FIBRO genomic record, causing the cell to differentiate into a 
FIBROBLAST CELL. 

Reference is now made to Fig. 5B, which is a simplified schematic diagram, 
illustrating cellular differentiation mediated by the € Genomic Records' concept. Fig. 5B depicts 
2 cells in an organism, designated CELL A and CELL B, each having a GENOME. It is 
appreciated that since CELL A and CELL B are cells in the same organism, the GENOME of 
CELL A is identical to that of CELL B. Despite having an identical GENOME, CELL A 
differentiates differently from CELL B, due to activation of different genomic records in these 
two cells. In CELL A the FIBRO GENOMIC RECORD is activated, causing CELL A to 
differentiate into a FIBROBLAST CELL, whereas in CELL B the BONE GENOMIC RECORD 
is activated, causing the CELL B to differentiate into a BONE CELL. The cause for activation of 
different genomic records in these two cells is the different maternal short RNA which they both 
received: CELL A received a maternal short RNA segment designated A' which activated 
genomic record FIBRO, whereas CELL B received a maternal short RNA segment designated B' 
which activated genomic record BONE. 
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Reference is now made to Fig. 6 which is a schematic diagram illustrating a 
'genomically p rogrammed c ell d ifferentiation' c oncept o f t he c onceptual m odel o f the p resent 
invention, addressing the genomic differentiation enigma. 

A cell designated CELL A divides into 2 cells designated CELL B and CELL C. 
CELL A, CELL B and CELL C each comprise a GENOME, which GENOME comprises a 
plurality of GENOMIC RECORDS. It is appreciated that since CELL A, CELL B and CELL C 
are cells in the same organism, the GENOME of these cells, and the GENOMIC RECORDS 
comprised therein, are identical 

As described above with reference to Fig. 5B, at its inception, CELL A receives a 
maternal short RNA segment, designated A' and marked by a broken line, which activates the 
FIBRO genomic record, thereby causing CELL A to differentiate into a FIBROBLAST CELL. 
However, Fig. 6 shows further details of the genomic records: each cell genomic record also 
comprises two short genomic sequences, referred to here as Daughter Cell Genomic Addresses. 
Blocks designated B and C are Daughter Cell Genomic Addresses of the FIBRO Genomic 
Record. At cell division, each parent cell transcribes two short RNA segments, corresponding to 
the two Daughter Cell Genomic Addresses of the Genomic Record of that parent cell, and 
transfers one to each of its two daughter cells. CELL A of Fig. 6 transcribes and transfers to its 
two respective daughter cells, two short RNA segments, outlined by a broken line and designated 
B' and C% corresponding to daughter cell genomic addresses designated B and C comprised in 
the FIBRO genomic record. 

CELL B therefore receives the above mentioned maternal short RNA segment 
designated B', which binds complementarity to genomic address designated B of genomic record 
BONE, thereby activating this genomic record, which in turn causes CELL B to differentiate 
into a BONE CELL. Similarly, CELL C receives the above mentioned maternal short RNA 
segment designated C\ which binds complementarity to genomic address designated C of 
genomic record CARTEL., thereby activating this genomic record, which in turn causes CELL C 
to differentiate into a CARTILAGE CELL. 

It is appreciated that the mechanism illustrated by Fig. 6 enables an unlimited 
lineage of cells to divide into daughter cells containing the same DNA, and to determine the cell- 
fate of these daughter cells. For example, when CELL B and CELL C divide into their respective 
daughter cells (not shown), they will transfer short RNA segments designated D' & E\ and F' & 
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G' respectively, to their respective daughter cells. The cell fate of each of these daughter cells 
would bed etermined b y t he i dentity o f t he m atemal s hort RNA se gment t hey receive, w hich 
would determine the genomic record activated. 

Reference is now made to Fig. 7 which is a schematic diagram illustrating a 
'genomically programmed cell-specific protein expression modulation' concept of the conceptual 
model of the present invention, addressing the genomic differentiation enigma. 

Cell A receives a maternal short RNA segment designated A', which activates a 
genomic record designated FIBRO, by anti-sense binding to a binding site 'header' of this 
genomic record, designated A. Genomic record FIBRO encodes 3 short RNA segments, 
designated 1, 2 and 4 respectively, which modulate expression of target genes designated 
GENE1, GENE2 and GENE4 respectively. Modulation of expression of these genes results in 
CELL A differentiating into a FIBROBLAST CELL. 

Reference is now made to Fig. 8, which is a simplified diagram describing each of 
a plurality of novel bioinformatically detected genes of the present invention, referred to here as 
Genomic Address Messenger (GAM) genes, which modulates expression of respective target 
genes thereof, the function and utility of which target genes is known in the art. 

GAM is a novel bioinformatically detected regulatory, non protein coding, micro 
RNA ( miRNA) gene. The m ethod b y w hich GAM i s d etected i s d escribed h ereinabove w ith 
reference to Figs. 8-15. 

GAM GENE and GAM TARGET GENE are human genes contained in the 
human genome. 

GAM GENE encodes a GAM PRECURSOR RNA. Similar to other miRNA 
genes, and unlike most ordinary genes, GAM PRECURSOR RNA does not encode a protein. 

GAM PRECURSOR RNA folds onto itself, forming GAM FOLDED 
PRECURSOR RNA, which has a two-dimensional 'hairpin structure'. As is well known in the 
art, this "hairpin structure', is typical of RNA encoded by miRNA genes, and is due to the fact 
that the nucleotide sequence of the first half of the RNA encoded by a miRNA gene is an 
accurately or partially inversed reversed sequence of the nucleotide sequence of the second half 
thereof. By inversed reversed is meant a sequence which is reversed and wherein each nucleotide 
is replaced by a complementary nucleotide, as is well known in the art (e.g. ATGGC is the 
complementary sequence of GCCAT). 
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An enzyme complex designated DICER COMPLEX, x dices^ the GAM FOLDED 
PRECURSOR RNA into GAM RNA, a single stranded -22 nt long RNA segment. As is known 
in the art, 'dicing' of a hairpin structured RNA precursor product into a short ~22nt RNA 
segment is catalyzed by an enzyme complex comprising an enzyme called Dicer together with 
other necessary proteins. 

TARGET GENE encodes a corresponding messenger RNA, GAM TARGET 
RNA. GAM TARGET RNA comprises three regions, as is typical of mRNA of a protein coding 
gene: a 5' untranslated region, a protein coding region and a 3' untranslated region, designated 
5'UTR, PROTEIN CODING and 3 X UTR respectively. 

GAM RNA binds complementarily to one or more target binding sites located in 
untranslated regions of GAM TARGET RNA. This complementary binding is due to the fact that 
the nucleotide sequence of GAM RNA is a partial or accurate inversed reversed sequence of the 
nucleotide sequence of each of the target binding sites. As an illustration, Fig. 8 shows three 
such target binding sites, designated BINDING SITE I, BINDING SITE II and BINDING SITE 
HI respectively. It is appreciated that the number of target binding sites shown in Fig. 8 is meant 
as an illustration only, and is not meant to be limiting GAM RNA may have a different number 
of target binding sites in untranslated regions of a GAM TARGET RNA. It is further appreciated 
that while Fig. 8 depicts target binding sites in the 3'UTR region, this is meant as an example 
only, these target binding sites may be located in the 3"UTR region, the 5^UTR region, or in both 
TXJTR and 5TOR regions. 

The complementary binding of GAM RNA to target binding sites on GAM 
TARGET RNA, such as BINDING SITE I, BINDING SITE H and BINDING SITE EI, inhibits 
translation of GAM TARGET RNA into GAM TARGET PROTEIN. GAM TARGET 
PROTEIN is therefore outlined by a broken line. 

It is appreciated that GAM TARGET GENE in fact represents a plurality of GAM 
target genes. The mRNA of each one of this plurality of GAM target genes comprises one or 
more target binding sites, each having a nucleotide sequence which is at least partly 
complementary to GAM RNA, and which when bound by GAM RNA causes inhibition of 
translation of respective one or more GAM target proteins. 

It is further appreciated by one skilled in the art that the mode of translational 
inhibition illustrated by Fig. 8 with specific reference to translational inhibition exerted by GAM 
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GENE on one or more TARGET GENE, is in fact common to other known miRNA genes, as is 
well known in the art. 

Nucleotide sequences of each of a plurality of GAM GENEs described by Fig. 8 
and their respective genomic source and chromosomal location are further described hereinbelow 
with reference to Table 1, hereby incorporated by reference. 

Nucleotide sequences of GAM PRECURSOR RNA, and a schematic 
representation of a predicted secondary folding of GAM FOLDED PRECURSOR RNA, of each 
of a plurality of GAM GENEs described by Fig. 8 are further described hereinbelow with 
reference to Table 2, hereby incorporated by reference. 

Nucleotide sequences of a Miced' GAM RNA of each of a plurality of GAM 
GENEs described by Fig. 8 are further described hereinbelow with reference to Table 3, hereby 
incorporated by reference. 

Nucleotide sequences of target binding sites, such as BINDING SITE-I, 
BINDING SITE-n and BINDING SITE-IE of Fig. 8, found on GAM TARGET RNA, of each of 
a plurality of GAM GENEs described by Fig. 8, and schematic representation of the 
complementarity of each of these target binding sites to each of a plurality of GAM RNA 
described by Fig. 8 are described hereinbelow with reference to Table 4, hereby incorporated by 
reference. 

It is appreciated that specific functions and accordingly utilities of each of a 
plurality of GAM GENEs described by Fig. 8 correlate with, and may be deduced from, the 
identity of the TARGET GENEs that each of said plurality of GAM GENEs binds and inhibits, 
and the function of each of said TARGET GENEs, as elaborated hereinbelow with reference to 
Table 5, hereby incorporated by reference. 

Studies establishing known functions of each of a plurality of TARGET GENEs of 
GAM GENEs of Fig. 8, and correlation of said each of a plurality of TARGET GENEs to known 
diseases are listed in Table 6, and are hereby incorporated by reference. 

The present invention discloses a novel group of genes, the GAM genes, 
belonging to the miRNA genes group, and for which a specific complementary binding has been 
determined. 
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Reference is now made to Fig. 9 which is a simplified block diagram illustrating a 
bioinformatic gene detection system capable of detecting genes of the novel group of genes of the 
present invention, which system is constructed and operative in accordance with a preferred 
embodiment of the present invention. 

A c enterpiece o f t he p resent i nvention i s a b ioinformatic gene d etection engine 
100, which is a preferred implementation of a mechanism capable of bioinformatically detecting 
genes of the novel group of genes of the present invention. 

The function of the bioinformatic gene detection engine 100 is as follows: it 
receives three types of input, expressed RNA data 102, sequenced DNA data 104, and protein 
function data 106, performs a complex process of analysis of this data as elaborated below, and 
based on this analysis produces output of a bioinformatically detected group of novel genes 
designated 108. 

Expressed RNA data 102 comprises published expressed sequence tags (EST) 
data, published mRNA data, as well as other sources of published RNA data. Sequenced DNA 
data 104 comprises alphanumeric data describing sequenced genomic data, which preferably 
includes annotation data such as location of known protein coding regions relative to the 
sequenced data. Protein function data 106 comprises scientific publications reporting studies 
which elucidated physiological function known proteins, and their connection, involvement and 
possible u tility i n t reatment a nd d iagnosis o f v arious d iseases. E xpressed RNA d ata 1 02 a nd 
sequenced DNA data 104 may preferably be obtained from data published by the National Center 
for Bioinformatics (NCBI) at the National Institute of Health (N1H) (Jenuth, 2000), as well as 
from various other published data sources. Protein function data 106 may preferably be obtained 
from any one of numerous relevant published data sources, such as the Online Mendelian 
Inherited D isease In M an ( OMIM(TM)) d atabase d eveloped b y Jo hn H opkins University, and 
also published by NCBI (2000). 

Prior to actual detection of bioinformatically detected novel genes 108 by the 
bioinformatic gene detection engine 100, a process of bioinformatic gene detection engine 
training & validation designated 110 takes place. This process uses the known miRNA genes as 
a training set (some 200 such genes have been found to date using biological laboratory means), 
to train the bioinformatic gene detection engine 100 to bioinformatically recognize miRNA-like 
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genes, and their respective potential target binding sites. Bioinformatic gene detection engine 
training & validation 1 10 is further described hereinbelow with reference to Fig. 10. 

The bioinformatic gene detection engine 100 comprises several modules which are 
preferably activated sequentially, and are described as follows: 

A non-coding genomic sequence detector 1 12 operative to bioinformatically detect 
non-protein coding genomic sequences. The non-coding genomic sequence detector 112 is 
further described herein below with reference to Figs. 1 1 A and 1 IB. 

A hairpin detector 114 operative to bioinformatically detect genomic 'hairpin- 
shaped x sequences, similar to GAM FOLDED PRECURSOR of Fig. 8. The hairpin detector 114 
is further described herein below with reference to Figs. 12A and 12B. 

A dicer-cut location detector 116 operative to bioinformatically detect the location 
on a hairpin shaped sequence which is enzymatically cut by DICER COMPLEX of Fig. 8. The 
dicer-cut location detector 1 16 is further described herein below with reference to Fig. 13 A. 

A target-gene binding-site detector 118 operative to bioinformatically detect target 
genes having binding sites, the nucleotide sequence of which is partially complementary to that 
of a given genomic s equence, s uch a s a s equence c ut b y DICER C OMPLEX o f F ig. 8 . T he 
target-gene binding-site detector 1 18 is further described hereinbelow with reference to Figs. 14A 
and 14B. 

A function & utility analyzer 120 operative to analyze function and utility of 
target genes, in order to identify target genes which have a significant clinical function and 
utility. The function & utility analyzer 120 is further described hereinbelow with reference to 
Fig. 15. 

Hardware implementation of the bioinformatic gene detection engine 100 is 
important, since significant computing power is preferably required in order to perform the 
computation of bioinformatic gene detection engine 100 in reasonable time and cost. For 
example, it is estimated that a using a powerful 8-processor server (e.g. DELL POWEREDGE 
(TM) 8450, 8 XEON (TM) 550MHz processors, 8 GB RAM), over 6 years (!) of computing time 
are required to detect all MIR genes in the human EST data, together with their respective 
binding sites. Various computer hardware and software configurations may be utilized in order 
to address this computation challenge, as is known in the art. A preferred embodiment of the 
present invention may preferably comprise a hardware configuration, comprising a cluster of one 
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hundred PCs (PENTIUM (TM) IV, 1.7GHz, with 40GB storage each), connected by Ethernet to 
12 servers (2-CPU, XEON (TM) 1.2-2.2GHz, with -200GB storage each), combined with an 8- 
processor server (8-CPU, Xeon 550Mhz w/ 8GB RAM) connected via 2 HBA fiber-channels to 
an EMC CLARHON (TM) 100-disks, 3.6 Terabyte storage device. A preferred embodiment of 
the present invention may also preferably comprise a software configuration which utilizes a 
commercial database software program, such as MICROSOFT (TM) SQL Server 2000. Using 
such preferred hardware and software configuration, may reduce computing time required to 
detect all MIR genes in the human EST data, and their respective binding sites, from 6 years to 
45 days. It is appreciated that the above mentioned hardware configuration is not meant to be 
limiting, and is given as an illustration only. The present invention may be implemented in a 
wide variety of hardware and software configurations. 

The p resent i nvention d iscloses 8 607 n ovel g enes o f t he G AM g roup o f g enes, 
which have been detected bioinformatically, as described hereinbelow with reference to Table 1 
through Table 6, and 1 096 novel genes of the GR group of genes, which have been detected 
bioinformatically, as described hereinbelow with reference to Table 7. Laboratory confirmation 
of 50 bioinformatically predicted genes of the GAM group of genes, and several 
bioinformatically predicted genes of the GR group of genes, is described hereinbelow with 
reference to Figs. 21 through 24C. 

Reference is now made to Fig. 10 which is a simplified flowchart illustrating 
operation of a mechanism for training a computer system to recognize the novel genes of the 
present invention. This mechanism is a preferred implementation of the bioinformatic gene 
detection engine training & validation 1 10 described hereinabove with reference to Fig. 9. 

BIOINFORMATIC GENE DETECTION ENGINE TRAINING & 
VALIDATION 110 of Fig. 9 begins by training the bioinformatic gene detection engine to 
recognize known miRNA genes, as designated by numeral 1 22. T his training step comprises 
HAIRPIN DETECTOR TRAINING & VALIDATION 124, further described hereinbelow with 
reference to Fig. 12A, DICER-CUT LOCATION DETECTOR TRAINING & VALIDATION 
126, further described hereinbelow with reference to Fig. 13A and 13B, and TARGET-GENE 
BINDING-SITE DETECTOR TRAINING & VALIDATION 128, further described hereinbelow 
with reference to Fig. 14 A. 



21 



WO 2004/044123 PCT/IL2003/000970 

Next, the BIOINFORMATIC GENE DETECTION ENGINE 100 is used to 
bioinformatically detect sample novel genes, as designated by numeral 130. Examples of sample 
novel genes thus detected are described hereinbelow with reference to Figs. 21 through 24C. 

Finally, wet lab experiments are preferably conducted in order to validate 
expression and preferably function of the sample novel genes detected by the BIOINFORMATIC 
GENE DETECTION ENGINE 100 in the previous step. An example of wet-lab validation of 
the above mentioned sample novel gene bioinformatically detected by the system is described 
hereinbelow with reference to Figs. 22A and 22B. 

Reference is now made to Fig. 1 1 A which is a simplified block diagram of a 
preferred implementation of the NON-CODING GENOMIC SEQUENCE DETECTOR 112 
described hereinabove with reference to Fig. 9. The NON-PROTEIN CODING GENOMIC 
SEQUENCE DETECTOR 112 of Fig. 9 preferably receives as input at least two types of 
published genomic data: EXPRESSED RNA DATA 102 and SEQUENCED DNA DATA 104. 
The EXPRESSED RNA DATA can include, among others, EST data, EST clusters data, EST 
genome alignment d ata and mRNA d ata. S ources for EXPRESSED RNA D ATA 1 02 include 
NCBI dbEST, NCBI UniGene clusters and mapping data, and TIGR gene indices. SEQUENCED 
DNA DATA 104 includes both sequence data (FASTA format files), and features annotation 
(GenBank file format) mainly from NCBI database. After its initial training, indicated by numeral 
134, and based on the above mentioned input data, the NON-PROTEIN CODING GENOMIC 
SEQUENCE DETECTOR 112 produces as output a plurality of NON-PROTEIN CODING 
GENOMIC SEQUENCES 136. Preferred operation of the NON-PROTEIN CODING 
GENOMIC SEQUENCE DETECTOR 112 is described hereinbelow with reference to Fig. 11B. 

Reference is now made to Fig. 11B which is a simplified flowchart illustrating a 
preferred operation of the NON-CODING GENOMIC SEQUENCE DETECTOR 112 of Fig. 9. 
Detection of NON-PROTEIN CODING GENOMIC SEQUENCES 136, generally preferably 
progresses in one of the following two paths: 

A first path for detecting NON-PROTEIN CODING GENOMIC SEQUENCES 
136 begins by receiving a plurality of known RNA sequences, such as EST data. Each RNA 
sequence is first compared to all known protein-coding sequences, in order to select only those 
RNA sequences which are non-protein coding, i.e. intergenic or intronic. This can preferably be 
performed by sequence comparison of the RNA sequence to known protein coding sequences, 
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using one of many alignment algorithms known in the art, such as BLAST. This sequence 
comparison to the DNA preferably also provides the localization of the RNA sequence on the 
DNA. 

Alternatively, selection of non-protein coding RNA sequences and their 
localization to the DNA can be performed by using publicly available EST clusters data and 
genomic mapping databases, such as UNIGENE database published by NCBI or TIGR database, 
in order map expressed RNA sequences to DNA sequences encoding them, to find the right 
orientation of EST sequences, and to exclude ESTs which map to protein coding DNA regions, 
as is well known in the art. Public databases, such as TIGR, may also be used to map an EST to a 
cluster of ESTs, assumed to be expressed as one piece, and is known in the art as Tentative 
Human Consensus. Publicly available genome annotation databases, such as NCBFs 
GENBANK, may also be used to deduce expressed intronic sequences. 

Optionally, an attempt may be made to 'expand' the non-protein RNA sequences 
thus found, by searching for transcription start and end signals, upstream and downstream of 
location of the RNA on the DNA respectively, as is well known in the art. 

A second path for detecting non-protein coding genomic sequences starts by 
receiving DNA sequences. The DNA sequences are parsed into non protein coding sequences, 
based on published DNA annotation data, by extracting those DNA sequences which are between 
known protein coding sequences. Next, transcription start and end signals are sought. If such 
signals are found, and depending on their 'strength', probable expressed non-protein coding 
genomic sequences are yielded. Such approach is especially useful for identifying novel GAM 
genes which are found in proximity to other known miRNA genes, or other wet-lab validated 
GAM genes. Since, as described hereinbelow with reference to Fig. 16, GAM genes are 
frequently found in clusters, therefore sequences near a known miRNA are more likely to contain 
novel GAM genes. Optionally, sequence orthology, i.e. sequences conservation in an 
evolutionary related species, may be used to select genomic sequences having higher probability 
of containing expressed novel GAM genes. 

Reference is now made to Fig. 12 A which is a simplified block diagram of a 
preferred implementation of the HAIRPIN DETECTOR 114 described hereinabove with 
reference to Fig. 9. 
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The goal of the HAIRPIN DETECTOR 1 14 is to detect 'haiipin' shaped genomic 
sequences, similar to those of known miRNA genes. As mentioned hereinabove with reference to 
Fig. 8, a 'haiipin' genomic sequence refers to a genomic sequence which 'folds onto itself 
forming a hairpin like shape, due to the fact that nucleotide sequence of the first half of the 
nucleotide sequence is an accurate or partial complementary sequence of the nucleotide sequence 
of its second half. 

The HAIRPIN DETECTOR 114 of Fig. 9 receives as input a plurality of NON- 
PROTEIN CODING GENOMIC SEQUENCES 136 of Fig. 11 A. After a phase of HAIRPIN 
DETECTOR TRAINING & VALIDATION 124 of Fig. 10, the HAIRPIN DETECTOR 114 is 
operative to detect and output 'hairpin shaped' sequences, which are found in the input NON- 
PROTEIN CODING GENOMIC SEQUENCES 138. The hairpin shaped sequences detected by 
the HAIRPIN DETECTOR 114 are designated HAIRPINS ON GENOMIC SEQUENCES 138. 
Preferred operation of the HAIPJTN DETECTOR 1 14 is described hereinbelow with reference to 
Fig. 12B. 

The phase of HAIRPE>J DETECTOR TRAINING & VALIDATION 124 is an 
iterative process of applying the HAIRPIN DETECTOR 1 14 to known haiipin shaped miRNA 
genes, calibrating the HAIRPIN DETECTOR 114 such that it identifies the training set of known 
hairpins, as well as sequences which are similar thereto. In a preferred embodiment of the 
present invention, THE HAIRPIN DETECTOR TRAINING & VALIDATION 124 trains and 
validates each of the steps of operation of the HAIRPIN DETECTOR 114, which steps are 
described hereinbelow with reference to Fig. 12B. 

The hairpin detector training and validation 124 preferably uses two sets of data: a 
training set of known miRNA genes, such as 440 miRNA genes of H. sapiens, M. musculus, C. 
elegans, C. Brigssae and D. Melanogaster, annotated in RFAM database (Griffiths- Jones 2003), 
and a large 'background set' of hairpins found in expressed non-protein coding genomic 
sequences, such as a set of 21,985 haiipins found in Tentative Human Concensus (THC) 
sequences in TIGR database. The 'background set' is expected to comprise some valid, 
previously undetected miRNA hairpins, and many haiipins which are not miRNA haiipins. 

In order to validate the performance of the HAIRPIN DETECTOR 114, a 
validation method is preferably used, which validation method is a variation on the k-fold cross 
validation method (Mitchell, 1997). This preferred validation method is devised to better cope 
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with the nature of the training set, which includes large families of similar and even identical 
miRNAs. The training set is preferably first divided into clusters of miRNAs such that any two 
miRNAs that belong to different clusters have an Edit Distance score (see Algorithms and 
Strings, Dan Gusfield, Cambridge University Press, 1997) of at least D=3, i.e. they differ by at 
least 3 editing operations. Next, the group of clusters is preferably divided into k sets. Then 
standard k-fold cross validation is preferably performed on this group of clusters, preferably 
using k=5, such that the members of each cluster are all in the training set or in the test set. It is 
appreciated that without the prior clustering, standard cross validation methods results in much 
higher performance of the predictors due to the redundancy of training examples, within the 
genome of a species and across genomes of different species. 

In a preferred embodiment of the present invention, using the abovementioned 
validation method, the efficacy of the HAIRPIN DETECTOR 114 is indeed validated: for 
example, when a similarity threshold is chosen such that 90% of the published miRNA-precursor 
hairpins are successfully predicted, only 7.6% of the 21,985 background hairpins are predicted to 
be miRNA-precursors, some of which may indeed be previously unknown miRNA precursors. 

Reference is now made to Fig. 12B which is a simplified flowchart illustrating a 
preferred operation of the HAlRPIN DETECTOR 1 14 of Fig. 9. 

A hairpin structure is a secondary structure, resulting from the nucleotide 
sequence pattern: the nucleotide sequence of the first half of the hairpin is a partial or accurate 
inversed reversed sequence of the nucleotide sequence of the second half thereof. Various 
methodologies are known in the art for prediction of secondary and tertiary hairpin structures, 
based on given nucleotide sequences. 

In a preferred embodiment of the present invention, the HAIRPIN DETECTOR 
114 initially calculates possible secondary structure folding patterns of a given one of the non- 
protein coding genomic sequences 136 and the respective energy of each of these possible 
secondary folding patterns, preferably using a secondary structure folding algorithm based on 
free-energy mi n imization, such as the MFOLD algorithm (Mathews et al., 1999), as is well 
known in the art. 

Next, the HAIRPIN DETECTOR 114 analyzes the results of the secondary 
structure folding, in order to determine the presence, and location of hairpin folding structures. A 
secondary structure folding algorithm, such as MFOLD algorithm, typically provides' as output a 
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listing of the base-pairing of the folded shape, i.e. a listing of each pair of connected nucleotides 
in the sequence. The goal of this second step is to asses this base-pairing listing, in order to 
determine if it describes a hairpin type bonding pattern. Preferably, each of the sequences that is 
determined to describe a hairpin structure is folded separately in order to determine its exact 
folding pattern and free-energy. 

The HAIRPIN DETECTOR 114 then assess those haiipin structures found by the 
previous step, comparing them to hairpins of known miRNA genes, using various characteristic 
hairpin features such as length of the hairpin and of its loop, free-energy and thermodynamic 
stability, amount and type of mismatched nucleotides, existence of sequence repeat-elements. 
Only hairpins that bear statistically significant resemblance to the training set of known miRNA 
hairpins, according to the abovementioned parameters are accepted. 

In a preferred embodiment of the present invention, similarity to the training set of 
known miRNA hairpins is determined using a 'similarity score' which is calculated using a 
weighted sum of terms, where each term is a function of one of the abovementioned hairpin 
features, and the parameters of each function are learned from the set of known hairpins, as 
described hereinabove with reference to hairpin detector training & validation 124. The weight 
of each term in the similarity score is optimized so as to achieve maximal separation between the 
distribution of similarity scores of hairpins which have been validated as miRNA-precursor 
hairpins, and the distribution of similarity scores of hairpins detected in the 'background set' 
mentioned hereinabove with reference to Fig. 12B, many of which are expected not to be 
miRNA-precursor hairpins. 

In another preferred embodiment of the present invention, the abovementioned 
DETERMINE IF SIMILAR TO KNOWN HAIRPIN-GENES step may may preferably be split 
into two stages. The first stage is a permissive filter that implements a simplified scoring method, 
based on a subset of the hairpin features described hereinabove, such as minimal length and 
maximal free energy. The second stage is more stringent, and a full calculation of the weighted 
sum of terms described hereinabove is performed. This second stage may preferably be 
performed only on the subset of hairpins that survived prior filtering stages of the hairpin- 
detector 114. 

Lastly, the HAIRPIN DETECTOR 114 attempts to select those hairpin structures 
which are as thermodynamically stable as the hairpins of known miRNA genes. This may 
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preferably be achieved in various manners. A preferred embodiment of the present invention 
utilizes the following methodology preferably comprising three logical steps: 

First, the HAIRPIN DETECTOR 114 attempts to group potential hairpins into 
'families' of closely related hairpins. As is known in the art, a free-energy calculation algorithm, 
typically provides multiple 'versions' each describing a different possible secondary structure 
folding pattern for the given genomic sequence, and the free energy of such possible folding. 
The HAIRPIN DETECTOR 114 therefore preferably assesses all hairpins found in each of the 
'versions', grouping hairpins which appear in different versions, but which share near identical 
locations into a common 'family' of hairpins. For example, all hairpins in different versions, the 
center of which hairpins is within 7 nucleotides of each other may preferably be grouped to a 
single 'family'. Hairpins may also be grouped to a single 'family' if the sequences of one or more 
hairpins are identical to, or are subsequences of, the sequence of another hairpin. 

Next, hairpin 'families' are assessed, in order to select only those families which 
represent hairpins that are as stable as those of known miRNA hairpins. Preferably only families 
which are represented in a majority of the secondary structure folding versions, such as at least in 
65% or 80% or 100% of the secondary structure folding versions, are considered stable. 

Finally, an attempt is made to select the most suitable hairpin from each selected 
family. For example, a hairpin which appears in more versions than other hairpins, and in 
versions the free-energy of which is lower, may be preferred. 

In another preferred embodiment of the present invention, hairpins with homology 
to other species, and clusters of thermodynamically stable hairpin are further favored. 

Reference is now made to Fig. 13A which is a simplified block diagram of a 
preferred implementation of the DICER-CUT LOCATION DETECTOR 116 described 
hereinabove with reference to Fig. 9. 

The goal of the DICER-CUT LOCATION DETECTOR 116 is to detect the 
location in which DICER COMPLEX of Fig. 8, comprising the enzyme Dicer, would 'dice' the 
given hairpin sequence, similar to GAM FOLDED PRECURSOR RNA, yielding GAM RNA 
both of Fig. 8. 

The DICER-CUT LOCATION DETECTOR 116 of Fig. 9 therefore receives as 
input a plurality of HAIRPEMS ON GENOMIC SEQUENCES 138 of Fig. 12A which were 
calculated by the previous step, and after a phase of DICER-CUT LOCATION DETECTOR 
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TRAINING & VALIDATION 126, is operative to detect a respective plurality of DICER-CUT 
SEQUENCES FROM HAIRPINS 140, one for each hairpin. 

In a preferred embodiment of the present invention, the DICER-CUT LOCATION 
DETECTOR 116 preferably uses standard machine learning techniques such as K nearest- 
neighbors, Bayesian networks and Support Vector Machines (SVM), trained on known dicer-cut 
locations of known miRNA genes in order to predict dicer-cut locations of novel GAM genes. 
The DICER-CUT LOCATION DETECTOR TRAINING & VALIDATION 126 is further 
described hereinbelow with reference to Fig. 13B. 

Reference is now made to Fig. 13B which is a simplified flowchart illustrating a 
preferred implementation of DICER-CUT LOCATION DETECTOR TRAINING & 
VALIDATION 126 of Fig. 10. 

The general goal of the DICER-CUT LOCATION D ETECTOR TRAINING & 
VALIDATION 126 is to analyze known hairpin shaped miRNA-precursors and their respective 
dicer-cut miRNA, in order to determine a common pattern to the dicer-cut location of known 
miRNA genes. Once such a common pattern is deduced, it may preferably be used by the 
DICER-CUT LOCATION DETECTOR 116, in detecting the predicted DICER-CUT 
SEQUENCES FROM HAIRPINS 140, from the respective HAIRPINS ON GENOMIC 
SEQUENCES 138, all of Fig. 13A. 

First, the dicer-cut location of all known miRNA genes is obtained and studied, so 
as to train the DICER-CUT LOCATION DETECTOR 116: for each of the known miRNA, the 
location of the miRNA relative to its hairpin-shaped miRNA-precursor is noted. 

The 5' and 3' ends of the dicer-cut location of each of the known miRNA genes is 
represented relative to the respective miRNA precursor hairpin, as well as to the nucleotides in 
each location along the hairpin. Frequency and identity of nucleotides and of nucleotide-pairing, 
and position of nucleotides and nucleotide pairing relative to the dicer-cut location in the known 
miRNA precursor hairpins is analyzed and modeled. In a preferred embodiment of the present 
invention, features learned from published miRNAs include: distance from hairpin's loop, 
nucleotide content, positional distribution of nucleotides and mismatched-nucleotides, and 
symmetry of mismatched-nucleotides. 

Different techniques are well known in the art of machine learning for analysis of 
existing pattern from a given 'framing set' of examples, which techniques are then capable, to a 
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certain degree, to detect similar patterns in other, previously unseen examples. Such machine 
learning techniques include, but are not limited to neural networks, Bayesian networks, Support 
Vector Machines (SVM), Genetic Algorithms, Markovian modeling, Maximum Liklyhood 
modeling, Nearest Neighbor algorithms, Decision trees and other techniques, as is well known in 
the art. 

The DICER-CUT LOCATION DETECTOR 116 preferably uses such standard 
machine learning techniques to predict either the 5' end or both the 5' and 3' ends of the miRNA 
excised, or 'diced' by the Dicer enzyme from the miRNA hairpin shaped precursor, based on 
known pairs of miRNA-precursors and their respective resulting miRNAs. The nucleotide 
sequences o f 4 40 p ublished m iRNA a nd t heir corresponding h airpin p recursors are p referably 
used for training and evaluation of the dicer-cut location detector module. 

Using the abovementioned training set, machine learning predictors, such as a 
Support Vector Machine (SVM) predictor, are implemented, which predictors test every possible 
nucleotide on a hairpin as a candidate for being the 5' end or the 3' end of a miRNA. Other 
machine learning predictors include predictors based on Nearest Neighbor, Bayesian modeling, 
and K-nearest-neighbor algorithms. The training set of the published miRNA precursor 
sequences is preferably used for training multiple separate classifiers or predictors, each of which 
produces a model for the 5' or 3' end of a miRNA relative to its hairpin precursor. The models 
take into account various miRNA properties such as the distance of the respective (3' or 5') end 
of the miRNA from the hairpin's loop, the nucleotides at its vicinity and the local 'bulge' (i.e. 
base-pair mismatch) structure. 

Performance of the resulting predictors, evaluated on the abovementioned 
validation set of 440 published miRNAs using k-fold cross validation (Mitchell, 1997) with k = 
3, is found to be as follows: in 70% of known miRNAs 5 '-end location is correctly determined by 
an SVM predictor within up to 2 nucleotides; a Nearest Neighbor (EDIT DISTANCE) predictor 
achieves 53% accuracy (233/440); a Two-Phased predictor that uses Baysian modeling (TWO 
PHASED) achieves 79% accuracy (348/440), when only the first phase is used, and 63% 
(277/440) when both phases are used; a K-nearest-neighbor predictor (FIRST-K) achieves 61% 
accuracy (268/440). The accuracies of all predictors are considerably higher on top scoring 
subsets of published miRNA. 
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Finally, in order to validate the efficacy and accuracy of the dicer-cut location 
detector 116, a sample of novel genes detected thereby is preferably selected, and validated by 
wet lab. Laboratory results validating the efficacy of the dicer-cut location detector 116 are 
described hereinbelow with reference to Figs. 21 through 24C. 



Reference i s n ow m ade to F ig. 1 3C w hich i s a s implified flowchart i Uustrating 
operation of DICER-CUT LOCATION DETECTOR 116 of Fig. 9, constructed and operative in 
accordance with a preferred embodiment of the present invention. 

The DICER CUT LOCATION DETECTOR 116 is a machine learning computer 
program module, which is trained on recognizing dicer-cut location of known miRNA genes, and 
based on this training, is operable to detect dicer cut location of novel GAM FOLDED 
PRECURSOR RNA. In a preferred embodiment of the present invention, the dicer-cut location 
module preferably utilizes machine learning algorithms, such as Support Vector Machine 
(SVM), Bayesian modeling, Nearest Neighbors, and K-nearest-neighbor, as is well known in the 
art. 

When assessing a novel GAM precursor, all 19-24 nucleotide long segments 
comprised in the GAM precursor are initially considered as 'potential GAMs', since the dicer-cut 
location is initially unknown. 

For each such potential GAM, its 5' end, or its 5' and 3' ends are scored by two or 
more recognition classifiers or predictors. 

In a preferred embodiment of the present invention, the DICER-CUT LOCATION 
DETECTOR 116 preferably uses a Support Vector Machine predictor trained on features such as 
distance from hairpin's loop, nucleotide content, positional distribution of nucleotides and 
mismatched-nucleotides, and symmetry of mismatched-nucleotides. 

In another preferred embodiment of the present invention, the DICER-CUT 
LOCATION DETECTOR 116 preferably uses an 'EDIT DISTANCE' predictor, which seeks 
sequences that are similar to those of published miRNAs, utilizing the Nearest Neighbor 
algorithm, where the similarity metric between two sequences is a variant of the edit distance 
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algorithm (Algorithms and Strings, Dan Gusfield, Cambridge University Press, 1997). This 
predictor is based on the observation that miRNAs tend to form clusters (Dostie, 2003), the 
members of which show marked sequence similarity to each other. 

In yet another preferred embodiment of the present invention, the DICER-CUT 
LOCATION DETECTOR 116 preferably uses a 'TWO PHASED' predictor, which predicts the 
dicer-cut location in two distinct phases: (a) selecting the double-stranded segment of the hairpin 
comprising the miRNA by naive Bayesian modeling (Mitchell, 1997), and (b) detecting which 
strand contains the miRNA by either naive or by K-nearest-neighbor modeling. The latter is a 
variant o f t he TJRST-K' p redictor d escribed h erein b elow, w ith p arameters o ptimized for t his 
specific task. The TWO PHASED' predictor may be operated in two modes: either utilizing only 
the first phase and thereby producing two alternative dicer-cut location predictions, or utilizing 
both phases and thereby producing only one final dicer-cut location. 

In still another preferred embodiment of the present invention, the DICER-CUT 
LOCATION DETECTOR 116 preferably uses a 'F1RST-K' predictor, which utilizes the K- 
nearest-neighbor algorithm. The similarity metric between any two sequences is 1- E/L , where L 
is a parameter, preferably 8-10 and E is the edit distance between the two sequences, taking into 
account only the first L nucleotides of each sequence. If the K-nearest-neighbor scores of two or 
more locations on the hairpin are not significantly different, these locations are further ranked by 
a Bayesian model, similar to the one described hereinabove. 

Scores of two or more of the abovementioned classifiers or predictors are 
integrated, yielding an integrated score for each 'potential GAM'. As an example, Fig. 13C 
illustrates integration of scores from two classifiers, a 3' end recognition classifier and a 5' end 
recognition classifier, the scores of which are integrated to yield an integrated score. In a 
preferred embodiment of the present invention, INTEGRATED SCORE of 13C preferably 
implements a 'best-of-breed' approach, accepting only 'potential GAMs' that score highly on one 
of the above mentioned EDIT DISTANCE', or 'TWO-PHASED'predictors. In this context, 
'high scores' means scores which have been demonstrated to have low false positive value when 
scoring known miRNAs. 

The INTEGRATED SCORE is then evaluated as follows: (a) the 'potential GAM' 
having the highest score is taken to be the most probable GAM, and (b) if the integrated score of 
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this 'potential GAM' is higher than a pre-defined threshold, then the potential GAM is accepted 
as the PREDICTED GAM. 

Reference is now made to Fig. 14A which is a simplified block diagram of a 
preferred implementation of the TARGET-GENE BINDING-SITE DETECTOR 118 described 
hereinabove with reference to Fig. 9. The goal of the TARGET-GENE BINDING-SITE 
DETECTOR 118 is to detect a BINDING SITE of Fig. 8, including binding sites located in 
untranslated regions of the RNA of a known gene, the nucleotide sequence of which BINDING 
SITE is a partial or accurate inversed reversed sequence to that of a GAM RNA of Fig. 8, thereby 
detennining that the above mentioned known gene is a target gene of GAM of Fig. 8. 

The TARGET-GENE BINDING-SITE DETECTOR 118 of Fig. 9 therefore 
receives as input a plurality of DICER-CUT SEQUENCES FROM HAIRPINS 140 of Fig. 13 A 
which were calculated by the previous step, and a plurality of POTENTIAL TARGET GENE 
SEQUENCES 142 which derive from SEQUENCED DNA DATA 104 of Fig. 9, and after a 
phase of TARGET-GENE BINDING-SITE DETECTOR TRAINING & VALIDATION 128 is 
operative to detect a plurality of POTENTIAL NOVEL TARGET-GENES HAVING BINDING 
SITE/S 144 the nucleotide sequence of which is a partial or accurate inversed reversed sequence 
to that of each of the plurality of DICER-CUT SEQUENCES FROM HAIRPINS 140. Preferred 
operation of the TARGET-GENE BINDING-SITE DETECTOR 118 is further described 
hereinbelow with reference to Fig. 14B. 

Reference is now made to Fig. 14B which is a simplified flowchart illustrating a 
preferred operation of the target-gene binding-site detector 1 18 of Fig. 9. 

In a preferred embodiment of the present invention, the target-gene binding-site 
detector 118 first uses a sequence comparison algorithm such as BLAST in order to compare the 
nucleotide sequence of each of the plurality of dicer-cut sequences from hairpins 140, to the 
potential target gene sequences 142, such a untranslated regions of known mRNAs, in order to 
find crude potential matches. Alternatively, the sequence comparison may preferably be 
performed using a sequence match search tool that is essentially a variant of the EDIT 
DISTANCE algorithm described hereinabove with reference to Fig. 13C, and the Nearest 
Neighbor algorithm (Mitchell, 1997). 

Results of the sequence comparison, performed by BLAST or other algorithms 
such as EDIT DISTANCE, are then filtered, preferably utilizing BLAST or EDIT DISTANCE 
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score, to results which are similar to those of known binding sites (e.g. binding sites of miRNA 
genes Lin-4 and Let-7 to target genes Lin-14, Lin-41, Lin 28 etc.). Next the binding site is 
expanded, checking if nucleotide sequenced immediately adjacent to the binding site found by 
the sequence comparison algorithm (e.g. BLAST or EDIT DISTANCE), may improve the match. 
Suitable binding sites, then are computed for free-energy and spatial structure. The results are 
analyzed, accepting only those binding sites, which have free-energy and spatial structure similar 
to that of known binding sites. Since known binding sites of known miRNA genes frequently 
have multiple adjacent binding sites on the same target RNA, accordingly binding sites which are 
clustered are strongly preferred. Binding sites found in evolutionarily conserved sequences may 
preferably also be preferred. 

For each candidate binding site a score, Binding Site Prediction Accuracy, is 
calculated which estimates their similarity of its binding to that of known binding sites. This 
score is based on GAM -binding site folding features including, but not limited to the free-energy 
, the total number and distribution of base pairs, the total number and distribution of unpaired 
nucleotides. 

In another preferred embodiment of the present invention binding sites are 
searched by a reversed process: s equences o f K(preferably 22)nucleotides o f the untranslated 
regions of the target gene are assessed as potential binding sites. A sequence comparison 
algorithm, such as BLAST or EDIT DISTANCE, is then used to search for partially or accurately 
complementary sequences elsewhere in the genome, which complementary sequences are found 
in known miRNA genes or computationally predicted GAM genes. Only complementary 
sequences, the complementarity of which withstands the spatial structure and free energy analysis 
requirements described above are accepted. Clustered binding sites are strongly favored, as are 
potential binding sites and potential GAM genes which occur in evolutionarily conserved 
genomic sequences. 

Target binding sites, identified by the TARGET-GENE BINDING-SITE 
DETECTOR 118, are divided into 3 groups: a) comprises binding sites that are exactly 
complementary to the predicted GAM. b) and c) comprise binding sites that are not exactly 
complementary to the predicted GAM: b) has binding sites with 0.9 < Binding Site Prediction 
Accuracy <= 1 and c) has binding sites with 0.8 < Binding Site Prediction Accuracy <= 0.9 . The 
average number of mismatching nucleotides in the alignment of predicted GAM and target 
binding site is smallest in category A and largest in category C. 
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In a preferred embodiment of the current invention a ranking of GAM to target 
gene binding is performed by calculating a score, Mir Target Accuracy. This score is the 
dominant group identifier of all binding sites of a specific GAM to a specific target gene UTR, 
where f a f dominates V and V dominates 'c\ 

In yet another preferred embodiment of the current invention a ranking of GAM to target 
gene binding is performed directly from the set of Binding Site Prediction Accuracies 
corresponding to all the binding sites of a specific GAM to a specific target gene UTR. This set 
of scores is sorted in descending order. The final score is a weighted sum of these scores where 
the weights are exponentially decreasing as a function of the rank. 

Reference is now made to Fig. 15 which is a simplified flowchart illustrating a 
preferred operation of the function & utility analyzer 120 described hereinabove with reference to 
Fig. 9. The goal of the function & utility analyzer 120 is to determine if a potential target gene is 
in fact a valid clinically useful target gene. Since a potential novel GAM gene binding a binding 
site in the UTR of a target gene is understood to inhibit expression of that target gene, and if that 
target gene is shown to have a valid clinical utility, then in such a case it follows that the potential 
novel gene itself also has a valid useful function - which is the opposite of that of the target gene. 

The function & utility analyzer 120 preferably receives as input a plurality of 
potential novel target genes having binding-site/s 144, generated by the target-gene binding-site 
detector 118, both of Fig. 14A. Each potential gene, is evaluated as follows: 

First, the system checks to see if the function of the potential target gene is 
scientifically well established. Preferably, this can be achieved bioinformatically by searching 
various published data sources presenting information on known function of proteins. Many such 
data sources exist and are published as is well known in the art. 

Next, for those target genes the function of which is scientifically known and is 
well documented, the system then checks if scientific research data exists which links them to 
known diseases. For example, a preferred embodiment of the present invention utilizes the 
OMIM(TM) database published byNCBI, which summarizes research publications relating to 
genes which have been shown to be associated with diseases. 
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Finally, the specific possible utility of the target gene is evaluated. W hile this 
process too may be facilitated by bioinformatic means, it might require manual evaluation of 
published scientific research regarding the target gene, in order to determine the utility of the 
target gene to the diagnosis and or treatment of specific disease. Only potential novel genes, the 
target-genes of which have passed all three examinations, are accepted as novel genes. 



Reference is now made to Fig. 16, which is a simplified diagram describing each 
of a plurality of novel bioinformatically detected regulatory genes, referred to here as Genomic 
Record (GR) genes, which encodes an 'operon-like' cluster of novel micro RNA-like genes, each 
of which in turn modulates expression of at least one target gene, the function and utility of 
which at least one target gene is known in the art. 

GR GENE is a novel bioinformatically detected regulatory, non protein coding, 
RNA gene. The method by which GR GENE was detected is described hereinabove with 
reference to Figs. 9-18, 

GR GENE encodes GR PRECURSOR RNA, an RNA molecule, typically several 
hundred nucleotides long. 

GR PRECURSOR RNA folds spatially, forming GR FOLDED PRECURSOR 
RNA. It is appreciated that GR FOLDED PRECURSOR RNA comprises a plurality of what is 
known in the art as 'hairpin' structures. These 'hairpin" structures are due to the fact that the 
nucleotide sequence of GR PRECURSOR RNA comprises a plurality of segments, the first half 
of each such segment having a nucleotide sequence which is at least a partial or accurate inversed 
reversed sequence of the second half thereof, as is well known in the art. 

GR FOLDED PRECURSOR RNA is naturally processed b y c ellular enzymatic 
activity into separate GAM precursor RNAs, herein schematically represented by GAM1 
FOLDED PRECURSOR RNA through GAM3 FOLDED PRECURSOR RNA, each of which 
GAM precursor RNAs being a hairpin shaped RNA segment, corresponding to GAM FOLDED 
PRECURSOR RNA of Fig.8. 

The above mentioned GAM precursor RNAs are diced by DICER COMPLEX of 
Fig.8, yielding respective short RNA segments of about 22 nucleotides in length, schematically 
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represented by GAM1 RNA through GAM3 RNA, each of which GAM RNAs corresponding to 
GAM RNA of Fig. 8. 

GAM1 RNA, GAM2 RNA and GAM3 RNA, each bind complementarity to binding sites located 
in untranslated r egions o f r espective t arget g enes, d esignated G AMI-TARGET RNA, G AM2- 
TARGET RNA and GAM3-TARGET RNA, respectively, which target binding site corresponds 
to a target binding site such as BINDING SITE I, BINDING SITE H or BINDING SITE m of 
Fig.8. This binding inhibits translation of the respective target proteins designated GAM1- 
TARGET PROTEIN, GAM2-TARGET PROTEIN and GAM3-TARGET PROTEIN 
respectively. 

It is appreciated that specific functions, and accordingly utilities, of each GR 
GENE of the present invention, correlates with, and may be deduced from, the identity of the 
target genes, which are inhibited by GAM RNAs comprised in the 'operon-like' cluster of said 
GR GENE, schematically represented by GAM1 TARGET PROTEIN through GAM3 TARGET 
PROTEIN. 

A listing of GAM GENEs comprised in each of a plurality of GR GENEs of Fig. 
16 is provided in Table 7, hereby incorporated by reference. Nucleotide sequences of each said 
GAM GENEs and their respective genomic source and chromosomal location are further 
described hereinbelow with reference to Table 1, hereby incorporated by reference. TARGET 
GENEs of each of said GAM GENEs are elaborated hereinbelow with reference to Table 4, 
hereby incorporated by reference. The functions of each of said TARGET GENEs and their 
association with various diseases, and accordingly the utilities of said each of GAM GENEs, and 
hence the functions and utilities of each of said GR GENEs of Fig. 16 is elaborated hereinbelow 
with reference to Table 5, hereby incorporated by reference. Studies establishing known 
functions of each of said TARGET GENEs, and correlation of each of said TARGET GENEs to 
known diseases are listed in table 6, and are hereby incorporated by reference. 

The present invention discloses 1096 novel genes of the GR group of genes, which 
have been detected bioinformatically, as elaborated hereinbelow with reference to Table 7. 
Laboratory confirmation of 2 genes of the GR group of genes is described hereinbelow with 
reference to Figs. 23A through 24C. 

In summary, the current invention discloses a very large number of novel GR 
genes, each of which encodes a plurality of GAM genes, which in turn may modulate expression 
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of a plurality of target proteins. It is appreciated therefore that the function of GR genes is in fact 
similar to that of the Genomic Records concept of the present invention addressing the 
differentiation enigma, described hereinabove with reference to Fig. 7. 

Reference is now made to Fig. 17 which is a simplified diagram illustrating a 
mode by which genes of a novel group of operon-like genes, described hereinabove with 
reference to Fig. 16 of the present invention, modulate expression of other such genes, in a 
cascading manner. 

GR1 GENE and GR2 GENE are two genes of the novel group of operon-like 
genes designated GR of Fig. 16. As is typical of genes of the GR group of genes, GR1 and GR2 
each encode a long RNA precursor, which in turn folds into a folded RNA precursor comprising 
multiple hairpin shapes, and is cut into respective separate hairpin shaped RNA segments, each of 
which RNA segments being diced to yield a gene of a group of genes designated GAM of Fig. 8. 
In this manner GR1 yields GAM1 RNA, GAM2 RNA and GAM3 RNA, and GR2 yields GAM4 
RNA, GAMS RNA and GAM6 RNA. 

As Fig. 17 shows, GAM3 RNA , which derives from GR1, binds a binding site 
located adjacent to GR2 GENE, thus modulating expression of GR2, thereby invoking expression 
of GAM4 RNA , GAM5 RNA and GAM6 RNA which derive from GR2. 

It is appreciated that the mode of modulation of expression presented by Fig. 17 
enables an unlimited "cascading effect" in which a GR gene comprises multiple GAM genes, each 
of which may modulate expression of other GR genes, each such GR gene comprising additional 
GAM genes, etc., whereby eventually certain GAM genes modulate expression of target proteins. 
This mechanism is in accord with the conceptual model of the present invention addressing the 
differentiation enigma, described hereinabove with specific reference to Figs. 6 and 7. 

Reference is now made to Fig. 18 which is a block diagram illustrating an 
overview of a methodology for finding novel genes and operon-like genes of the present 
invention, and their respective functions. 

According to a preferred embodiment of the present invention, the methodology to 
finding novel genes of the present invention and their function comprises of the following major 
steps: 

First, genes of the novel group of genes of the present invention, referred to here 
as GAM genes, are located and their function elicited by detecting target proteins they bind and 
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the function of those target proteins, as described hereinabove with reference to Figs. 9 through 
15. 

Next, genes of a novel group of operon-like genes of the present invention, 
referred to here as GR genes, are located, by locating clusters of proximally located GAM genes, 
based on the previous step. 

Consequently, the hierarchy of GR and GAM genes is elicited: binding sites for 
non-protein-binding GAM genes comprised in each GR gene found are sought adjacent to other 
GR genes. When found, such a binding site indicates that the connection between the GAM and 
the GR the expression of which it modulates, and thus the hierarchy of the GR genes and the 
GAM genes they comprise. 

Lastly, the function of GR genes and GAM genes which are "high" in the 
hierarchy, i .e. G AM genes w hich m odulate expression o f o ther GR g enes r ather t han d irectly 
modulating expression of target proteins, may be deduced. A preferred approach is as follows: 
The function of protein-modulating GAM genes is deducible from the proteins which they 
modulate, provided that the function of these target proteins is known. The function of "higher' 
GAM genes may be deduced by .comparing the function of protein-modulating GAM genes, with 
the hierarchical relationships by which the 'higher" GAM genes are connected to the protein- 
modulating GAM genes. For example, given a group of several protein-modulating GAM genes, 
which collectively cause a protein expression pattern typical of a certain cell-type, then a 'higher" 
GAM gene is sought which modulates expression of GR genes which perhaps modulate 
expression of other genes which eventually modulate expression of the given group of protein- 
modulating GAM genes. The "higher" GAM gene found in this manner is taken to be responsible 
for differentiation of that cell-type, as per the conceptual model of the invention described 
hereinabove with reference to Fig. 6. 

Reference is now made to Fig. 19 which is a block diagram illustrating different 
utilities of genes of the novel group of genes of the present invention referred to here as GAM 
genes and GR genes. 

The present invention discloses a first plurality of novel genes referred to here as 
GAM genes, and a second plurality of operon-like genes referred to here as GR genes, each of the 
GR genes encoding a plurality of GAM genes. The present invention further discloses a very 
large number of known target-genes, which are bound by, and the expression of which is 
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modulated by each of the novel genes of the present invention. Published scientific data 
referenced by the present invention provides specific, substantial, and credible evidence that the 
above mentioned target genes modulated by novel genes of the present invention, are associated 
with various diseases. Specific novel genes of the present invention, target genes thereof and 
diseases associated therewith, are described hereinbelow with reference to Tables 1 through 7. It 
is therefore appreciated that a function of GAM genes and GR genes of the present invention is 
modulation of expression of target genes related to known diseases, and that therefore utilities of 
novel genes of the present invention include diagnosis and treatment of the above mentioned 
diseases. Fig. 19 describes various types of diagnostic and therapeutic utilities of novel genes of 
the present invention. 

A utility of novel genes of the present invention is detection of GAM genes and of 
GR genes. It is appreciated that since GAM genes and GR genes modulate expression of disease 
related target genes, that detection of expression of GAM genes in clinical scenarios associated 
with said diseases is a specific, substantial and credible utility. Diagnosis of novel genes of the 
present invention may preferably be implemented by RNA expression detection techniques, 
including but not limited to biochips, as is well known in the art. Diagnosis of expression of 
genes of the present invention may be useful for research purposes, in order to further understand 
the connection between the novel genes of the present invention and the above mentioned related 
diseases, for disease diagnosis and prevention purposes, and for monitoring disease progress. 

Another utility of novel genes of the present invention is anti-GAM gene therapy, 
a mode of therapy which allows up regulation of a disease related target-gene of a novel GAM 
gene of the present invention, by lowering levels of the novel GAM gene which naturally inhibits 
expression of that target gene. This mode of therapy is particularly useful with respect to target 
genes which have been shown to be under-expressed in association with a specific disease. Anti- 
GAM gene therapy is further discussed hereinbelow with reference to Figs. 20A and 20B. 

A further utility of novel genes of the present invention is GAM replacement 
therapy, a mode of therapy which achieves down regulation of a disease related target-gene of a 
novel GAM gene of the present invention, by raising levels of the GAM gene which naturally 
inhibits expression of that target gene. This mode of therapy is particularly useful with respect to 
target genes which have been shown to be over-expressed in association with a specific disease. 
GAM replacement therapy involves introduction of supplementary GAM gene products into a 
cell, or stimulation of a cell to produce excess GAM gene products. GAM replacement therapy 
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may preferably be achieved b y transfecting c ells with an artificial DNA molecule encoding a 
GAM gene, which causes the cells to produce the GAM gene product, as is well known in the art. 
Yet a further utility of novel genes of the present invention is modified GAM therapy. Disease 
conditions are likely to exist, in which a mutation in a binding site of a GAM geneprevents 
natural GAM gene to effectively bind inhibit a disease related target-gene, causing up regulation 
of that target gene, and thereby contributing to the disease pathology. In such conditions, a 
modified GAM gene is designed which effectively binds the mutated GAM binding site, i.e. is an 
effective anti-sense of the mutated GAM binding site, and is introduced in disease effected cells. 
Modified GAM therapy is preferably achieved by transfecting cells with an artificial DNA 
molecule encoding the modified GAM gene, which causes the cells to produce the modified 
GAM gene product, as is well known in the art. 

An additional utility of novel genes of the present invention is induced cellular 
differentiation therapy. As aspect of the present invention is finding genes which determine 
cellular differentiation, as described hereinabove with reference to Fig. 18. Induced cellular 
differentiation therapy comprises transfection of cell with such GAM genes thereby determining 
their differentiation as desired. It is appreciated that this approach may be widely applicable, 
inter aha as a means for auto transplantation harvesting cells of one cell-type from a patient, 
modifying their differentiation as desired, and then transplanting them back into the patient. It is 
further appreciated that this approach may also be utilized to modify cell differentiation in vivo, 
by transfecting cells in a genetically diseased tissue with a cell-differentiation determining GAM 
gene, thus stimulating these cells to differentiate appropriately. 

Reference is now made to Figs. 20A and 20B, simplified diagrams which when 
taken together illustrate anti-GAM gene therapy m entioned hereinabove with reference to Fig. 
19. A utility of novel genes of the present invention is anti-GAM gene therapy, a mode of 
therapy which allows up regulation of a disease related target-gene of a novel GAM gene of the 
present invention, by lowering levels of the novel GAM gene which naturally inhibits expression 
of that target gene. Fig. 20A shows a normal GAM gene, inhibiting translation of a target gene of 
GAM gene, by binding to a BINDING SITE found in an untranslated region of GAM TARGET 
RNA, as described hereinabove with reference to Fig. 8. 

Fig. 20B shows an example of anti-GAM gene therapy. ANTI-GAM RNA is 
short artificial RNA molecule the sequence of which is an anti-sense of GAM RNA. Anti-GAM 
treatment comprises transfecting diseased cells with ANTI-GAM RNA, or with a DNA encoding 
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thereof. The ANTI-GAM RNA binds the natural GAM RNA, thereby preventing binding of 
natural GAM RNA to its BINDING SITE. This prevents natural translation inhibition of GAM 
TARGET RNA by GAM RNA, thereby up regulating expression of GAM TARGET PROTEIN. 

It is appreciated that anti-GAM gene therapy is particularly useful with respect to 
target genes which have been shown to be under-expressed in association with a specific disease. 

Furthermore, anti-GAM therapy is particularly useful, since it may be used in 
situations in which technologies known in the art as RNAi and siRNA can not be utilized. As in 
known in the art, RNAi and siRNA are technologies which offer means for artificially inhibiting 
expression of a target protein, by artificially designed short RNA segments which bind 
complementarily to mRNA of said target protein. However, RNAi and siRNA can not be used to 
directly upregulate translation of target proteins. 

Reference is now made to Fig. 21, which is a table summarizing laboratory 
validation results that validate efficacy of the bioinformatic gene detection engine 100 of Fig. 9. 
In order to assess efficacy of the bioinformatic gene detection engine 100, novel genes predicted 
thereby are preferably divided into 4 DETECTION ACCURACY GROUPS (first column), 
designated A through D, ranking GAMS from the most probable GAMs to the least probable 
GAMs, using the scores of HAIRPIN DETECTOR 114 and DICER-CUT LOCATION 
DETECTOR 1 16 as follows: 

Group A: The score of the HAIRPIN-DETECTOR is above 0.7, the overall score 
of the two-phased predictor is above 0.55, and the score of the second phase of the two-phased 
predictor is above 0.75. Group B: The score of the EDIT-DISTANCE predictor is equal or above 
17. lii both groups A and B one miRNA is predicted for each hairpin. Group C: The score of the 
HAIRPIN-DETECTOR is between 0.6 and 0.7, and the overall score of the two-phased predictor 
is above 0.55. Group D: The score of the HAIRPIN-DETECTOR is between 0.3 and 0.6, and the 
overall score of the two-phased predictor is above 0.55. In both groups C and D, for each hairpin 
both sides of the double stranded window are given as output, and are examined in the lab. The 
groups are mutually exclusive, i.e. in groups A, C and D all hairpins score less than 17 in the 
EDIT-DISTANCE predictor. It is appreciated that the division into groups is not exhaustive: 344 
of the 440 published hairpins, and 4747 of the 8607 novel GAMs of the present invention, belong 
to one of the groups. 
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Sample novel bioinfoimatically predicted genes, of each of these groups are sent 
to the laboratory for validation (third column), and the number (fourth column) and percent (fifth 
column) of successful validation of predicted genes is noted for each of the groups, as well as 
overall (bottom line). The number of novel genes explicitly specified by present invention 
belonging to each of the four groups is noted (sixth column), as is the number of novel genes 
found in intergenic and intronic RNA (seventh and eighth columns). Expected positives (eighth 
column) is deduced by multiplying the % success in the lab validation by the number of novel 
GAMs detected in intronic and intergenic RNA databases. 

It is appreciated that the present invention comprises 4747 novel GAM genes, 
which fall into one of these four detection accuracy groups, and that the bioinformatic prediction 
is substantiated by a group of 50 novel GAM genes validated by laboratory means, out of 147 
genes which were tested in the lab, resulting in validation of an overall 34% accuracy. The top 
group demonstrated 50% accuracy. Pictures of test-results of specific genes in the 
abovementioned four groups, as well as the methodology used for validating the expression of 
predicted genes is elaborated hereinbelow with reference to Fig. 22. 

It is further appreciated that failure to detect a gene in the lab does not necessarily 
indicate a mistaken bioinformatic prediction. Rather, it may be due to technical sensitivity 
limitation of the lab test, or because the gene is not expressed in the tissue examined, or at the 
development phase tested. 

It is still further appreciated that in general these findings are in agreement with 
the expected bioinformatic accuracy, as describe hereinabove with reference to Fig. 13B: 
assuming 80% accuracy of the hairpin detector 114 and 80% accuracy of the dicer-cut location 
detector 1 16 and 80% accuracy of the lab validation, this would result in 50% overall accuracy of 
the genes validated in the lab. 

Reference is now made to Fig. 22 which is a picture of laboratory results 
validating the expression of 25 novel genes detected by the bioinformatic gene detection engine 
100, in the four detection accuracy groups A through D described hereinabove with reference to 
Fig. 21. 

Each row in Fig. 22, designated A through D, correlates to a corresponding one of 
the four detection accuracy groups A-D, described hereinabove with reference to Fig. 21. In each 
row, pictures of several genes validated by hybridization of PCR-product southern-blots, are 

42 



WO 2004/044123 PCT7IL2003/000970 

provided, each corresponding to a specific GAM gene, as elaborated hereinbelow. These PCR- 
product hybridization pictures are designated 1 through 5 in the A group, 1 through 12 in the B 
group, 1 through 9 in the C group, and 1 through 4 in the D group. In each PCR hybridization 
picture, 2 lanes are seen: the test lane, designated and the control lane, designated For 
convenience of viewing the results, all PCR-product hybridization pictures of Fig. 22 have been 
shrunk x4 vertically. It is appreciated that for each of the tested genes, a clear hybridization band 
appears in the test ('+') lane, but not in the control ('-') lane. 

Specifically, Fig 22 shows pictures of PCR-product hybridization validation by 
southern-blot, the methodology of which is described hereinbelow, to the following novel GAM 
genes: 

DETECTION ACCURACY GROUP A: (1) GAM8297.1; (2) GAM5346.1; 
(3) GAM281.1; (4) GAM8554.1; and (5)GAM2071.1. 

DETECTION ACCURACY GROUP B: (1) GAM7553.1; (2) GAM5385.1; 
(3) GAM5227.1; (4) GAM7809.1; (5) GAM 1032.1; (6) GAM3431.1; (7) GAM7933.1; (8) 
GAM3298.1.;(9) GAM116.1; (10) GAM3418.1 (later published by other researchers as MIR23); 
(11) GAM3499.1; (12) GAM3027.1; (13) GAM7080.1; (14) GAM895.1; and (15) GAM2608.1, 
(16) GAM20, and (17) GAM21. 

DETECTION ACCURACY GROUP C: (1) GAM3770.1; (2) GAM1338.1; 
(3)GAM7957.1;(4)GAM391.1; (5) GAM 8678.1; (6) GAM2033.1; (7)GAM7776.1; (8) 
GAM8145.1; and 

(9) GAM 633.1. 

DETECTION ACCURACY GROUP D: (1) GAM19; (2) GAM8358.1; 
(3) GAM3229.1; and (4)GAM7052.1. 

In addition to the PCR detection, the following GAMs were cloned and 
sequenced: GAM1338.1, GAM7809.1, GAM116.1, GAM3418.1 (later published by other 
researchers as MIR23), GAM3499.1, GAM3027.1, GAM7080.1, and GAM21. 

The PCR-product hybridization validation methodology used is briefly described 
as follows. In order to validate the expression of predicted novel GAM genes, and assuming that 
these novel genes are probably expressed at low concentrations, a PCR product cloning approach 
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was set up through the following strategy: two types of cDNA libraries designated "One tailed" 
and "Ligation" were prepared from frozen HeLa S100 extract (4c Biotech, Belgium) size 
fractionated RNA. Essentially, Total S100 RNA was prepared through an SDS-Proteinase K 
incubation followed by an acid Phenol-Chloroform purification and Isopropanol precipitation. 
Alternatively, total HeLa RNA was also used as starting material for these libraries. 

Fractionation was done by loading up to 500ug per YM100 Amicon Microcon 
column (Millipore) followed by a 500g centrifugation for 40 minutes at 4»C. Flowthrough 
"YM100" RNA consisting of about % of the total RNA was used for library preparation or 
fractionated further by loading onto a YM30 Amicon Microcon column (Millipore) followed by 
a 13,500g centrifugation for 25 minutes at 4°C. Flowthrough "YM30" was used for library 
preparation as is and consists of less than 0.5% of total RNA. For the both the "ligation" and the 
"One-tailed" libraries RNA was dephosphorilated and ligated to an RNA (lowercase)-DNA 
(UPPERCASE) hybrid 5'-phosphorilated, 3' idT blocked 3'-adapter (5'-P- 
uuuAACCGCATTCTC-idT-3' Dharmacon # P-002045-01-05) (as elaborated in Elbashir et al 
2001) resulting in ligation only of RNase m type cleavage products. 3'-Ligated RNA was 
excised and purified from a half 6%, half 13% polyacrylamide gel to remove excess adapter with 
a Nanosep 0.2uM centrifugal device (Pall) according to instructions, and precipitated with 
glycogen and 3 volumes of Ethanol. Pellet was resuspended in a ininimal volume of water. 
For the "ligation" library a DNA (UPPERCASE)-RNA (lowercase) hybrid 5'-adapter (5'- 
TACTAATACGACTCACTaaa-3 ' Dharmacon # P-002046-01-05) was ligated to the 3'-adapted 
RNA, reverse transcribed with "EcoRI-RT" : (5'. 

GACTAGCTGGAATTCAAGGATGCGGTTAAA-3 '), PCR amplified with two external 
primers essentially as in Elbashir et al 2001 except that primers were "EcoRI-RT" and 'TstI 
Fwd" (5 ' -CAGCC AACGCTGCAGATACGACTC ACTAAA-3 '). This PCR product was used as 
a template for a second round of PCR with one hemispecific and one external primer or with two 
hemispecific primers. 

For the "One tailed" library the 3 '-Adapted RNA was annealed to 20pmol primer 
"EcoRI RT" by heating to 70°C and cooling 0.1»C/sec to 30°C and then reverse transcribed with 
Superscript H RT (According to instructions, Invitrogen) in a 20ul volume for 10 alternating 5 
minute cycles of 37°C and 45°C. Subsequently, RNA was digested with lul 2M NaOH, 2mM 
EDTA at 65°C for 10 minutes. cDNA was loaded on a polyacrylamide gel, excised and gel- 
purified from excess primer as above (invisible, judged by primer run alongside) and 
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resusperided in 13ul of water. Purified cDNA was then oligo-dC tailed with 400U of recombinant 
terminal transferase (Roche molecular biochemicals), lul lOOuM dCTP, lul 15mM CoCl 2 , and 
4nl reaction buffer, to a final volume of 20ul for 15 minutes at 37°C. Reaction was stopped with 
2ul 0.2M EDTA and 15ul 3M NaOAc pH 5.2. Volume was adjusted to 150ul with water, Phenol 
: Bromochloropropane 10:1 extracted and subsequently precipitated with glycogen and 3 
volumes of Ethanol. C-tailed cDNA was used as a template for PCR with the external primers 
'T3-PstBsg(G/I) 18 " " (5 ,_ 

AATTAACCCTCACTAAAGGCTGCAGGTGCAGGIGGGIIGGGIIGGGIIGN-3' where I 
stands for Inosine and N for any of the 4 possible deoxynucleotides), and with "EcoRI Nested" 
(5 '-GGAATTCAAGGATGCGGTTA-3 '). This PCR product was used as a template for a second 
round of PCR with one hemispecific and one external primer or with two hemispecific primers. 

Hemispecific primers were constructed for each predicted GAM by an in-house 
program designed to choose about half of the 5' or 3' sequence of the GAM corresponding to a 
TM° of about 30°-34°C constrained by an optimized 3' clamp, appended to the cloning adapter 
sequence (for "One-tailed" libraries 5 '-GGNNGGGNNG on the 5' end of the GAM, or 
TTTAACCGCATC-3 ' on the 3' end of the GAM. For "Ligation" libraries the same 3' adapter 
and 5 '-CGACTCACTAAA on the 5' end). Consequently, a fully complementary primer of a 
TM° higher than 60°C was created covering only one half of the GAM sequence permitting the 
unbiased elucidation by sequencing of the other half. 

CONFIRMATION OF GAM SEQUENCE AUTHENTICITY OF PCR PRODUCTS: 

SOUTHERN BLOT: PCR-product sequences were confirmed by southern blot 
(Southern EM. Biotechnology.l992;24: 122-39.(1975)) and hybridization with DNA 
oligonucleotide probes synthesized against predicted GAMs. Gels were transferred onto a 
Biodyne PLUS 0.45pm, (Pall) positively charged nylon membrane and UV cross-linked. 
Hybridization was performed overnight with DIG-labeled probes at 42°C in DIG Easy-Hyb 
buffer (Roche). Membranes were washed twice with 2xSSC and 0.1% SDS for 10 min. at 42°C 
and then washed twice with 0.5xSSC and 0.1% SDS for 5 min at 42°C. The membrane was then 
developed by using a DIG luminescent detection kit (Roche) using anti-DIG and CSPD reaction, 
according to the manufacturer's protocol. All probes were prepared according to the 
manufacturers (Roche Molecular Biochemicals) protocols: Digoxigenin (DIG) labeled antisense 
transcripts was prepared from purified PCR products using a DIG RNA labeling kit with T3 
RNA polymerase. DIG labeled PCR was prepared by using a DIG PCR labeling kit. 3'-DIG- 
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tailed oligo ssDNA antisense probes, containing DIG-dUTP and dATP at an average tail length 
of 50 nucleotides were prepared from lOOpmole oligonucleotides with the DIG Oligonucleotide 
Labeling Kit. 

CLONING: PCR products were inserted into pGEM-T (Promega) or pTZ57 
(MBI Fermentas), transformed into competent JM109 E. coli (Promega) and sown on LB-Amp 
plates with IPTG/Xgal. White and light-blue colonies were transferred to duplicate gridded 
plates, one of which was blotted onto a membrane (Biodyne Plus, Pall) for hybridization with 
DIG tailed oligo probes (according to instructions, Roche) corresponding to the expected GAM. 
Plasmid DNA from positive colonies was sequenced. 

Reference is now made to Fig. 23 A, which is a schematic representation of a novel 
human GR gene, herein designated GR12731, located on chromosome 9, comprising 2 known 
MIR genes - M1R24 and MIR23, and 2 novel GAM genes, herein designated GAM22 and 
GAM1 16, all marked by solid black boxes. Fig. 23A also schematically illustrates 6 non-GAM 
hairpin sequences, and one non-hairpin sequence, all marked by white boxes, and serving as 
negative controls. By 'non-GAM hairpin sequences' is meant sequences of a similar length to 
known MIR PRECURSOR sequences, which form hairpin secondary folding pattern similar to 
MIR PRECURSOR hairpins, and yet which are assessed by the bioinformatic gene detection 
engine 100 not to be valid GAM PRECURSOR hairpins. It is appreciated that Fig. 23A is a 
simplified schematic representation, reflecting only the order in which the segments of interest 
appear relative to one another, and not a proportional distance between the segments. 

Reference is now made to Fig. 23B, which is a schematic representation of 
secondary folding of each of the MIRs and GAMs of GR GR12731 - MIR24, MIR23, GAM22 
and GAM116, and of the negative control non-GAM hairpins, herein designated N2, N3, N116, 
N4, N6 and N7. NO is a non-hairpin control, of a similar length to that of known MIR 
PRECURSOR hairpins. It is appreciated that the negative controls are situated adjacent to and in 
between real MIR and GAM genes, and demonstrates similar secondary folding patterns to that 
of known MIRs and GAMs. 

Reference is now made to Fig. 23C, which is a picture of laboratory results of a 
PCR test upon a YM100 "ligation'Mibrary, utilizing specific primer sets directly inside the 
boundaries of the hairpins. Due to the nature of the library the only PCR amplifiable products 
can result from RNaselH type enzyme cleaved RNA, as expected for legitimate hairpin 
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precursors presumed to be produced by DROSHA (Lee et al, Nature 425 415-419, 2003). Fig 
23C demonstrates expression of hairpin precursors of known MIR genes - MIR23 and MIR24, 
and of novel bioinfonnatically detected GAM22 and GAM116 genes predicted bioinformatically 
by a system constructed and operative in accordance with a preferred embodiment of the present 
invention. Fig. 23C also shows that none of the 7 controls (6 hairpins designated N2, N3, N23, 
N4, N6 and N7 and 1 non-hairpin sequence designated NO) were expressed. Nl 16 is a negative 
control sequence partially overlapping GAM1 16. 

In the picture, test lanes including template are designated and the control lane is designated 
It is appreciated that for each of the tested hairpins, a clear PCR band appears in the test ('+*) 
lane, but not in the control ('-') lane. 

Figs. 23A through 23C, when taken together validate the efficacy of the 
bioinformatic gene detection engine in: (a) detecting known MIR genes; (b) detecting novel 
GAM genes which are found adjacent to these MIR genes, and which despite exhaustive prior 
biological e fforts a nd b ioinformatic d etection e fforts, w ent undetected; ( c) d iscerning b etween 
GAM (or MIR) PRECURSOR hairpins, and non-GAM hairpins. 

It is appreciated that the ability to discern GAM-hairpins from non-GAM-hairpins 
is very significant in detecting GAM genes, since hairpins in general are highly abundant in the 
genome. Other MIR prediction programs have not been able to address this challenge 
successfully. 

Reference is now made to Fig. 24A which is an annotated sequence of an EST 
comprising a novel gene detected by the gene detection system of the present invention. Fig. 
24A shows the nucleotide sequence of a known human non-protein coding EST (Expressed 
Sequence Tag), identified as EST72223. The EST72223 clone obtained from TIGR database 
(Kirkness and Kerlavage, 1997) was sequenced to yield the above 705bp transcript with a 
polyadenyl tail. It is appreciated that the sequence of this EST comprises sequences of one known 
miRNA gene, identified as M1R98, and of one novel GAM gene, referred to here as GAM25, 
detected by the bioinformatic gene detection system of the present invention and described 
hereinabove with reference to Fig. 9. 
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The sequences of the precursors of the known MIR98 and of the predicted 
GAM25 are in bold, the sequences of the established miRNA 98 and of the predicted miRNA 
GAM25 are underlined. 

Reference is now made to Figs. 24B and 24C that are pictures of laboratory 
results, which when taken together demonstrate laboratory confirmation of expression of the 
bioinformatically detected novel gene of Fig. 24 A. 

I In two parallel experiments, an enzymatically synthesized capped, EST72223 

RNA transcript, was incubated with Hela SI 00 lysate for 0 minutes, 4 hours and 24 hours. RNA 
was subsequently harvested, run on a denaturing polyacrylamide gel, and reacted with a 102nt 
and a 145nt antisense MIR98 and GAM25 transcript probes respectively. The Northern blot 
results of these experiments demonstrated processing of EST72223 RNA by Hela lysate (lanes 2- 
4, in 24B and 24C), into ~80bp and ~22bp segments, which reacted with the MIR98 probe (24B), 
and into ~100bp and ~24bp segments, which reacted with the GAM25 probe (24C). These results 
demonstrate the processing of EST72223 by Hela lysate into MIR98 and GAM25. It is also 
appreciated from Fig. 24C (lane 1) that Hela lysate itself reacted with the GAM25 probe, in a 
number of bands, including a ~100bp band, indicating that GAM25-precursor is endogenously 
expressed in Hela cells. The presence of additional bands, higher than lOObp in lanes 5-9 
probably corresponds to the presence of nucleotide sequences in Hela lysate, which contain the 
GAM25 sequence. 

In addition, in order to demonstrate the kinetics and specificity of the processing 
of MIR98 and GAM25 miRNA precursors into their respective miRNA's, transcripts of MIR98 
and of the bioinformatically predicted GAM25, were similarly incubated with Hela SI 00 lysate, 
for 0 minutes, 30 minutes, 1 hour and 24 hours, and for 24 hours with the addition of EDTA, 
added to inhibit Dicer activity, following which RNA was harvested, run on a polyacrylamide gel 
and reacted with MIR98 and GAM25 probes. Capped transcripts were prepared for in-vitro RNA 
cleavage assays with T7 RNA polymerase including a m 7 G(5 f )ppp(5 f )G-capping reaction using 
the m Message m Machine k it ( Ambion). P urified P CR p roducts w ere u sed a s t emplate for t he 
reaction. These were amplified for each assay with specific primers containing a T7 promoter at 
the 5* end and a T3 RNA polymerase promoter at the 3'end. Capped RNA transcripts were 
incubated at 30°C in supplemented, dialysis concentrated, Hela SI 00 cytoplasmic extract (4C 
Biotech, Seneffe, Belgium). The Hela S100 was supplemented by dialysis to a final concentration 
of 20mM Hepes, lOOmM KC1, 2.5mM MgCl 2 , 0.5mM DTT, 20% glycerol and protease inhibitor 
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cocktail tablets (Complete mini Roche Molecular Biochemicals). After addition of all 
components, final concentrations were lOOmM capped target RNA, 2mM ATP, 0.2mM GTP, 
500U/ml RNasin, 25ug/ml creatine kinase, 25mM creatine phosphate, 2.5mM DTT and 50% 
S100 extract. Proteinase K, used to enhance Dicer activity (Zhang H, Kolb FA Brondani V, Billy 
E, Filipowicz W. Human Dicer preferentially cleaves dsRNAs at their termini without a 
requirement for ATP.EMBO J. 2002 Nov l;21(21):5875-85) was dissolved in 50mM Tris-HCl 
pH 8, 5mM CaCh, and 50% glycerol, was added to a final concentration of 0.6 mg/ml. Cleavage 
reactions were stopped by the addition of 8 volumes of proteinase K buffer (200Mm Tris-Hcl, pH 
7.5, 25m M EDTA 300mM NaCl, and 2% SDS) and incubated at 65°C for 15min at different 
time points (0, 0.5, 1, 4, 24h) and subjected to phenol/chloroform extraction. Pellets were 
dissolved in water and kept frozen. Samples were analyzed on a segmented half 6%, half 13% 
polyacrylamide 1XTBE-7M Urea gel. 

The Northern blot results of these experiments demonstrated an accumulation of a 
~22bp segment which reacted with the MIR98 probe, and of a ~24bp segment which reacted with 
the GAM25 probe, over time (lanes 5-8). Absence of these segments when incubated with EDTA 
(lane 9), which is known to inhibit Dicer enzyme (Zhang et al., 2002), supports the notion that 
the processing of MER98 and GAM25 miRNAs from their precursors is mediated by Dicer 
enzyme, found in Hela lysate. The molecular sizes of EST72223, MIR-98 and GAM25 are 
indicated by arrows. 

Transcript products were 705nt (EST72223), 102nt (MIR98), 125nt (GAM25) 
long.EST72223 was PCR amplified with T7-EST 72223 forward primer: 

5 '-TAATACGACTCACTATAGGCCCTTATTAGAGGATTCTGCT-3 ' 
and T3-EST72223 reverse primer: 



MIR98 was PCR amplified using EST72223 as a template with T7MIR98 forward primer: 

5-TAATACGACTCACTATAGGGTGAGGTAGTAAGTTGTATTGTT-3' 

and T3MIR98 reverse primer: 
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5 , -AATTAACCCTCACTAAAGGGAAAGTAGTAAGTTGTATAGTT-3\ 



GAM25 was PCR amplified using EST72223 as a template with GAM25 forward primer: 5'- 
GAGGCAGGAGAATTGCTTGA- 3' and T3-EST72223 reverse primer: 

5 5 - AATTAACCCTCACTAAAGGCCTGAGAC AGAGTCTTGCTC-3 * . 

To validate the identity of the band shown by the lower arrow in fig. 24C, a RNA band 
parallel to a marker of 24 base was excised from the gel and cloned as in Elbashir et al (2001) 
and sequenced. 90 clones corresponded to the sequence of GAM25, three corresponded to 
GAM25* (the opposite arm of the hairpin with a 1-3 nucleotide 3' overhang) and two to the 
hairpin-loop. 

GAM25 was also validated endogenously by sequencing from both sides from a HeLa 
YM100 total-RNA "ligation" libraries, utilizing hemispecific primers as detailed in fig22. 

Taken together, these results validate the presence and processing of a novel MIR 
gene product, GAM25, which was predicted bioinformatically. The processing of this novel gene 
product, by Hela lysate .from EST72223, through its precursor, to its final form was similar to 
that observed for known gene, MIR98. 

It is appreciated that the data presented in Figs. 24A, 24B and 24C when taken 
together validate the function of the bioinformatic gene detection engine 100 of Fig. 9. Fig. 24A 
shows a novel GAM gene bioinformatically detected by the bioinformatic gene detection engine 
100, and Figs. 24B and 24C show laboratory confirmation of the expression of this novel gene. 
This is in accord with the engine training and validation methodology described hereinabove with 
reference to Fig. 10. 

It is appreciated by persons skilled in the art that the present invention is not 
limited by what has been particularly shown and described hereinabove. Rather the scope of the 
present invention includes both combinations and subcombinations of the various features 
described hereinabove as well as variations and modifications which would occur to persons 
skilled in the art upon reading the specifications and which are not in the prior art. 
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DETAILED DESCRIPTION OF LARGE TABLES 

Table 1 comprises data relating to the source and location of novel GAM genes of 
the present invention, and contains the following fields: 

GENE NAME Rosetta Genomics Ltd. gene nomenclature (see below) 

GAM SEQ-ID GAM Seq-ID, as in the Sequence Listing 

PRECUR SEQ-ID GAM precursor Seq-ID, as in the Sequence Listing 

ORGANISM Abbreviated (hsa = Homo sapiens) 

CHR Chromosome encoding the GAM gene 

CHROMOSOME OFFSET Offset of GAM precursor sequence on chromosome 
SOURCEREF-ID Accession number of source sequence 

SOURCE_OFFSET Offset of GAM precursor sequence on source sequence 

SRC Source-type of GAM precursor sequence (see below) 

GAM ACC GAM Prediction Accuracy Group (see below); 



Table 2 comprises data relating to GAM precursors of novel GAM genes of the 
present invention, and contains the following fields: 

GENE NAME Rosetta Genomics Ltd. gene nomenclature (see below) 

PRECUR SEQ-ID GAM precursor Seq-ID, as in the Sequence Listing 

PRECURSOR SEQUENCE GAM precursor nucleotide sequence (5' to 3 X ) 

FOLDED-PRECURSOR Schematic representation of the GAM folded precursor, 

beginning 5' end (beginning of upper row) to 3" end 
(beginning of lower row), where the hairpin loop is 
positioned at the right part of the draw. 

SRC Source-type of GAM precursor sequence (see below) 

GAM ACC GAM Prediction Accuracy Group (see below); 
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Table 3 comprises data relating to GAM genes of the present invention, and 
contains the following fields: 

Table3: Genes 



GENE NAME 
GAM SEQ-ID; 
GENE_SEQUENCE 
PRECUR SEQ-ID 
SOURCEJUEF-ID 
SRC 

GAM ACC 



Rosetta Genomics Ltd. gene nomenclature (see below) 
GAM Seq-ID, as in the Sequence Listing 
Sequence (5' to 3 X ) of the mature, Miced' GAM gene 
GAM precursor Seq-ID, as in the Sequence Listing 
Accession number of the source sequence 
Source-type of GAM precursor sequence (see below) 
GAM Prediction Accuracy Group (see below); 



Table 4 comprises data relating to target-genens and binding sites of GAM genes 
of the present invention, and contains the following fields: 



GENE NAME 
GAM SEQ-ID; 
TARGET 
#BS 

TARGET SEQ-ID 
TARGET REF-ID 
UTR 

UTR OFFSET 



Rosetta Genomics Ltd. gene nomenclature (see below) 

GAM Seq-ID, as in the Sequence Listing 

GAM target protein name 

Number of binding sites of GAM onto Target 

Target binding site Seq-ID, as in the Sequence Listing 

Target accession number (GenBank) 

Untranslated region of binding site/s (3* or 5') 

Offset of GAM binding site relative to UTR 

TAR-BINDING-SITE-SEQ Nucleotide sequence (5' to 3') of the target binding site 

BINDING-SITE-DRAW Schematic representation of the binding site, upper row 

present 5" to 3" sequence of the GAM, lower row present 3' 
to 5" sequence of the target. 

SRC Source-type of GAM precursor sequence (see below) 
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GAM ACC 
BS ACC 
TAR ACC 
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GAM Prediction Accuracy Group (see below); 

Binding-Site Accuracy Group (see below) 

Target Accuracy Group (see below); 



Table 5 comprises data relating to functions and utilities of novel GAM genes of 
the present invention, and contains the following fields: 



GENE NAME 
TARGET 

GENE_SEQUENCE 

GENE-FUNCTION 

SRC 

GAM ACC 
TAR ACC 
TAR DIS 



Rosetta Genomics Ltd. gene nomenclature (see below) 
GAM target protein name 

Sequence (5' to 3') of the mature, x diced % GAM gene 
Description of the GAM functions and utilities 
Source-type of GAM precursor sequence (see below) 
GAM Prediction Accuracy Group (see below) 
Target Accuracy Group (see below) 
Target Accuracy Group (see below); 



Table 6 comprises a bibliography of references supporting the functions and 
utilities of novel GAM genes of the present invention, and contains the following fields: 

Table6: Gene Function References- Bibliography. 



GENE NAME 
TARGET 
REFERENCES 
SRC 

GAM ACC 
TAR ACC 



Rosetta Genomics Ltd. gene nomenclature (see below) 
GAM target protein name 
list of references relating to the target gene, 
Source-type of GAM precursor sequence (see below) 
GAM Prediction Accuracy Group (see below) 
Target Accuracy Group (see below); and 
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Table 7 comprises data relating to novel GR genes of the present invention, and 
contains the following fields: 

GENE NAME Rosetta Genomics Ltd. GR gene nomenclature 

SOURCE START, OFFSET Start-offset of GR gene relative to source sequence 

SOURCE ENDJ3FFSET End-offset of GR gene relative to source sequence 

SOURCE_REF-ID Accession number of the source sequence 

GAMS JD'S_IN_GR List of the GAM genes in the GR cluster 

SRC Source-type of GAM precursor sequence (see below) 

GR ACC GR Prediction Accuracy Group (see below). 

The following conventions and abbreviations are used in the tables: 

GENE NAME is a RosettaGenomics Ltd. gene nomenclature. All GAMs are 
designated by GAMx.l or GAMx.2 where x is the unique SEQ-ID. If the GAM precursor has a 
single prediction for GAM, it is designated by GAMx. 1 . Otherwise, the higher accuracy GAM 
prediction is designated by GAMx.l and the second is designated by GAMx.2. 

SRC is a field indicating the type of source in which novel genes were detected, as 
one of the following options: (1) TIGR Intergenic, (3) EST or Unigene Intron Intergenic, (4) 
TIGR Intron, (6) DNA Intergenic, (7) DNA Intron, (8) DNA Exon. Sequences are based on 
NCBI Build33 of the human genome. TIGR source is based on "Tentative Human Consensus" 
(THC) The Institute for Genomic Research which are not found in mRNA Intron/Exon according 
to NCBI GenBank genome annotation. 

GAM ACC (GAM Prediction Accuracy Group) of gene prediction system: A- 
very high accuracy, B- high accuracy,C- moderate accuracy, D-low accuracy, as described 
hereinbelow with reference to Fig.21. 

BS ACC (Binding-Site Accuracy Group) indicates accuracy of total GAM-target 
binding prediction, considering the number of binding sites a GAM has on the target's UTR; A- 
very high accuracy, B- high accuracy, C- moderate accuracy, as described hereinbelow 

with reference to Fig.l4B. 
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TAR ACC (Target Accuracy Group) indicates accuracy of target binding site 
prediction, A- very high accuracy, B- high accuracy, C- moderate accuracy, as described 
hereinbelow with reference to Fig.l4B. 

TAR DIS (Target Disease Relation Group) 'A 1 indicates if the target gene is 
known to have a specific causative relation to a specific known disease, based on the OMIM 
database. It is appreciated that this is a partial classification emphasizing genes which are 
associated with * single gene* diseases etc. All genes of the present invention ARE associated 
with various diseases, although not all are in c A' status. 

GR ACC (GR Prediction Accuracy Group) indicates the maximum gene 
prediction accuracy among GAM genes of the cluster, A- very high accuracy, B- high accuracy, 
C- moderate accuracy, as described hereinbelow with reference to Fig.l4B. 
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